####################################################
# Steven R. Gehrke
# Northern Arizona University
# 22-0338: Bicycling Accessibility (CRANC 2.0)
# Task: Plots (Study 1: Access)
# Created: 07.13.2023
# Updated: 07.25.2023
####################################################

####################
### DATA SET-UP
####################

# Load libraries
# install.packages(c("tidyverse","foreign"))
packages <- c("tidyverse","foreign")
lapply(packages, require, character.only=TRUE); rm(packages)

# Set directory
dat_dir <- "F:/22-0338_BicyclingAccessibility/"
dat_dir2 <- "H:/_projects/22-0338_BicyclingAccessibility/"

# Import accessibility data sets
dat_ibc_e <- read_csv(paste(dat_dir, "dat_access_ibc_e.csv", sep="_data/_tabular/_outputs/"))
dat_eac_e <- read_csv(paste(dat_dir, "dat_access_eac_e.csv", sep="_data/_tabular/_outputs/"))
dat_saf_e <- read_csv(paste(dat_dir, "dat_access_saf_e.csv", sep="_data/_tabular/_outputs/"))
dat_ibc_m <- read_csv(paste(dat_dir, "dat_access_ibc_m.csv", sep="_data/_tabular/_outputs/"))
dat_eac_m <- read_csv(paste(dat_dir, "dat_access_eac_m.csv", sep="_data/_tabular/_outputs/"))
dat_saf_m <- read_csv(paste(dat_dir, "dat_access_saf_m.csv", sep="_data/_tabular/_outputs/"))
dat_ibc_s <- read_csv(paste(dat_dir, "dat_access_ibc_s.csv", sep="_data/_tabular/_outputs/"))
dat_eac_s <- read_csv(paste(dat_dir, "dat_access_eac_s.csv", sep="_data/_tabular/_outputs/"))
dat_saf_s <- read_csv(paste(dat_dir, "dat_access_saf_s.csv", sep="_data/_tabular/_outputs/"))

# Color palette
brew <- c("#e41a1c","#377eb8","#4daf4a","#984ea3","#ff7f00","#ffff33","#a65628","#f781bf")

####################
### EMPLOYMENT
####################

# Summarize employment access by cyclist type
ibc_e <- dat_ibc_e %>% group_by(mpo) %>%
  summarise(e_05=mean(s_00_05,na.rm=T),e_10=mean(s_05_10,na.rm=T),e_15=mean(s_10_15,na.rm=T),e_20=mean(s_15_20,na.rm=T),
            e_25=mean(s_20_25,na.rm=T),e_30=mean(s_25_30,na.rm=T),e_35=mean(s_30_35,na.rm=T),e_40=mean(s_35_40,na.rm=T),
            e_45=mean(s_40_45,na.rm=T),e_50=mean(s_45_50,na.rm=T),e_55=mean(s_50_55,na.rm=T),e_60=mean(s_55_60,na.rm=T),
            e_05_lts=mean(s_00_05_lt,na.rm=T),e_10_lts=mean(s_05_10_lt,na.rm=T),e_15_lts=mean(s_10_15_lt,na.rm=T),e_20_lts=mean(s_15_20_lt,na.rm=T),
            e_25_lts=mean(s_20_25_lt,na.rm=T),e_30_lts=mean(s_25_30_lt,na.rm=T),e_35_lts=mean(s_30_35_lt,na.rm=T),e_40_lts=mean(s_35_40_lt,na.rm=T),
            e_45_lts=mean(s_40_45_lt,na.rm=T),e_50_lts=mean(s_45_50_lt,na.rm=T),e_55_lts=mean(s_50_55_lt,na.rm=T),e_60_lts=mean(s_55_60_lt,na.rm=T))
eac_e <- dat_eac_e %>% group_by(mpo) %>%
  summarise(e_05=mean(s_00_05,na.rm=T),e_10=mean(s_05_10,na.rm=T),e_15=mean(s_10_15,na.rm=T),e_20=mean(s_15_20,na.rm=T),
            e_25=mean(s_20_25,na.rm=T),e_30=mean(s_25_30,na.rm=T),e_35=mean(s_30_35,na.rm=T),e_40=mean(s_35_40,na.rm=T),
            e_45=mean(s_40_45,na.rm=T),e_50=mean(s_45_50,na.rm=T),e_55=mean(s_50_55,na.rm=T),e_60=mean(s_55_60,na.rm=T),
            e_05_lts=mean(s_00_05_lt,na.rm=T),e_10_lts=mean(s_05_10_lt,na.rm=T),e_15_lts=mean(s_10_15_lt,na.rm=T),e_20_lts=mean(s_15_20_lt,na.rm=T),
            e_25_lts=mean(s_20_25_lt,na.rm=T),e_30_lts=mean(s_25_30_lt,na.rm=T),e_35_lts=mean(s_30_35_lt,na.rm=T),e_40_lts=mean(s_35_40_lt,na.rm=T),
            e_45_lts=mean(s_40_45_lt,na.rm=T),e_50_lts=mean(s_45_50_lt,na.rm=T),e_55_lts=mean(s_50_55_lt,na.rm=T),e_60_lts=mean(s_55_60_lt,na.rm=T))
saf_e <- dat_saf_e %>% group_by(mpo) %>%
  summarise(e_05=mean(s_00_05,na.rm=T),e_10=mean(s_05_10,na.rm=T),e_15=mean(s_10_15,na.rm=T),e_20=mean(s_15_20,na.rm=T),
            e_25=mean(s_20_25,na.rm=T),e_30=mean(s_25_30,na.rm=T),e_35=mean(s_30_35,na.rm=T),e_40=mean(s_35_40,na.rm=T),
            e_45=mean(s_40_45,na.rm=T),e_50=mean(s_45_50,na.rm=T),e_55=mean(s_50_55,na.rm=T),e_60=mean(s_55_60,na.rm=T),
            e_05_lts=mean(s_00_05_lt,na.rm=T),e_10_lts=mean(s_05_10_lt,na.rm=T),e_15_lts=mean(s_10_15_lt,na.rm=T),e_20_lts=mean(s_15_20_lt,na.rm=T),
            e_25_lts=mean(s_20_25_lt,na.rm=T),e_30_lts=mean(s_25_30_lt,na.rm=T),e_35_lts=mean(s_30_35_lt,na.rm=T),e_40_lts=mean(s_35_40_lt,na.rm=T),
            e_45_lts=mean(s_40_45_lt,na.rm=T),e_50_lts=mean(s_45_50_lt,na.rm=T),e_55_lts=mean(s_50_55_lt,na.rm=T),e_60_lts=mean(s_55_60_lt,na.rm=T))

### IBC Cyclist Type

# Graphic set-up
png(file=paste(dat_dir2, "plot_ibc-employment.png", sep="/_analysis/_plots/"), bg="white", width=9, height=9, units='in', res=750)
plot(c(0,60), c(0,50000), type="n", xlab="Travel time (minutes)", ylab="Job accessibility", xaxt="n", cex=1.2)
# Plot points and lines
# CYMPO
lines(5, ibc_e$e_05_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(10, ibc_e$e_10_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(15, ibc_e$e_15_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(20, ibc_e$e_20_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(25, ibc_e$e_25_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(30, ibc_e$e_30_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(35, ibc_e$e_35_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(40, ibc_e$e_40_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(45, ibc_e$e_45_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(50, ibc_e$e_50_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(55, ibc_e$e_55_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(60, ibc_e$e_60_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
segments(0, 0, 5, ibc_e$e_05_lts[1], col=brew[1], lwd=1.5)
segments(5,  ibc_e$e_05_lts[1], 10, ibc_e$e_10_lts[1], col=brew[1], lwd=1.5)
segments(10, ibc_e$e_10_lts[1], 15, ibc_e$e_15_lts[1], col=brew[1], lwd=1.5)
segments(15, ibc_e$e_15_lts[1], 20, ibc_e$e_20_lts[1], col=brew[1], lwd=1.5)
segments(20, ibc_e$e_20_lts[1], 25, ibc_e$e_25_lts[1], col=brew[1], lwd=1.5)
segments(25, ibc_e$e_25_lts[1], 30, ibc_e$e_30_lts[1], col=brew[1], lwd=1.5)
segments(30, ibc_e$e_30_lts[1], 35, ibc_e$e_35_lts[1], col=brew[1], lwd=1.5)
segments(35, ibc_e$e_35_lts[1], 40, ibc_e$e_40_lts[1], col=brew[1], lwd=1.5)
segments(40, ibc_e$e_40_lts[1], 45, ibc_e$e_45_lts[1], col=brew[1], lwd=1.5)
segments(45, ibc_e$e_45_lts[1], 50, ibc_e$e_50_lts[1], col=brew[1], lwd=1.5)
segments(50, ibc_e$e_50_lts[1], 55, ibc_e$e_55_lts[1], col=brew[1], lwd=1.5)
segments(55, ibc_e$e_55_lts[1], 60, ibc_e$e_60_lts[1], col=brew[1], lwd=1.5)
#FMPO
lines(5, ibc_e$e_05_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(10, ibc_e$e_10_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(15, ibc_e$e_15_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(20, ibc_e$e_20_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(25, ibc_e$e_25_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(30, ibc_e$e_30_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(35, ibc_e$e_35_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(40, ibc_e$e_40_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(45, ibc_e$e_45_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(50, ibc_e$e_50_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(55, ibc_e$e_55_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(60, ibc_e$e_60_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
segments(0, 0, 5, ibc_e$e_05_lts[2], col=brew[2], lwd=1.5)
segments(5,  ibc_e$e_05_lts[2], 10, ibc_e$e_10_lts[2], col=brew[2], lwd=1.5)
segments(10, ibc_e$e_10_lts[2], 15, ibc_e$e_15_lts[2], col=brew[2], lwd=1.5)
segments(15, ibc_e$e_15_lts[2], 20, ibc_e$e_20_lts[2], col=brew[2], lwd=1.5)
segments(20, ibc_e$e_20_lts[2], 25, ibc_e$e_25_lts[2], col=brew[2], lwd=1.5)
segments(25, ibc_e$e_25_lts[2], 30, ibc_e$e_30_lts[2], col=brew[2], lwd=1.5)
segments(30, ibc_e$e_30_lts[2], 35, ibc_e$e_35_lts[2], col=brew[2], lwd=1.5)
segments(35, ibc_e$e_35_lts[2], 40, ibc_e$e_40_lts[2], col=brew[2], lwd=1.5)
segments(40, ibc_e$e_40_lts[2], 45, ibc_e$e_45_lts[2], col=brew[2], lwd=1.5)
segments(45, ibc_e$e_45_lts[2], 50, ibc_e$e_50_lts[2], col=brew[2], lwd=1.5)
segments(50, ibc_e$e_50_lts[2], 55, ibc_e$e_55_lts[2], col=brew[2], lwd=1.5)
segments(55, ibc_e$e_55_lts[2], 60, ibc_e$e_60_lts[2], col=brew[2], lwd=1.5)
#LHMPO
lines(5, ibc_e$e_05_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(10, ibc_e$e_10_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(15, ibc_e$e_15_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(20, ibc_e$e_20_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(25, ibc_e$e_25_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(30, ibc_e$e_30_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(35, ibc_e$e_35_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(40, ibc_e$e_40_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(45, ibc_e$e_45_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(50, ibc_e$e_50_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(55, ibc_e$e_55_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(60, ibc_e$e_60_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
segments(0, 0, 5, ibc_e$e_05_lts[3], col=brew[3], lwd=1.5)
segments(5,  ibc_e$e_05_lts[3], 10, ibc_e$e_10_lts[3], col=brew[3], lwd=1.5)
segments(10, ibc_e$e_10_lts[3], 15, ibc_e$e_15_lts[3], col=brew[3], lwd=1.5)
segments(15, ibc_e$e_15_lts[3], 20, ibc_e$e_20_lts[3], col=brew[3], lwd=1.5)
segments(20, ibc_e$e_20_lts[3], 25, ibc_e$e_25_lts[3], col=brew[3], lwd=1.5)
segments(25, ibc_e$e_25_lts[3], 30, ibc_e$e_30_lts[3], col=brew[3], lwd=1.5)
segments(30, ibc_e$e_30_lts[3], 35, ibc_e$e_35_lts[3], col=brew[3], lwd=1.5)
segments(35, ibc_e$e_35_lts[3], 40, ibc_e$e_40_lts[3], col=brew[3], lwd=1.5)
segments(40, ibc_e$e_40_lts[3], 45, ibc_e$e_45_lts[3], col=brew[3], lwd=1.5)
segments(45, ibc_e$e_45_lts[3], 50, ibc_e$e_50_lts[3], col=brew[3], lwd=1.5)
segments(50, ibc_e$e_50_lts[3], 55, ibc_e$e_55_lts[3], col=brew[3], lwd=1.5)
segments(55, ibc_e$e_55_lts[3], 60, ibc_e$e_60_lts[3], col=brew[3], lwd=1.5)
#MAG
lines(5, ibc_e$e_05_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(10, ibc_e$e_10_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(15, ibc_e$e_15_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(20, ibc_e$e_20_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(25, ibc_e$e_25_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(30, ibc_e$e_30_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(35, ibc_e$e_35_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(40, ibc_e$e_40_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(45, ibc_e$e_45_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(50, ibc_e$e_50_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(55, ibc_e$e_55_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(60, ibc_e$e_60_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
segments(0, 0, 5, ibc_e$e_05_lts[4], col=brew[4], lwd=1.5)
segments(5,  ibc_e$e_05_lts[4], 10, ibc_e$e_10_lts[4], col=brew[4], lwd=1.5)
segments(10, ibc_e$e_10_lts[4], 15, ibc_e$e_15_lts[4], col=brew[4], lwd=1.5)
segments(15, ibc_e$e_15_lts[4], 20, ibc_e$e_20_lts[4], col=brew[4], lwd=1.5)
segments(20, ibc_e$e_20_lts[4], 25, ibc_e$e_25_lts[4], col=brew[4], lwd=1.5)
segments(25, ibc_e$e_25_lts[4], 30, ibc_e$e_30_lts[4], col=brew[4], lwd=1.5)
segments(30, ibc_e$e_30_lts[4], 35, ibc_e$e_35_lts[4], col=brew[4], lwd=1.5)
segments(35, ibc_e$e_35_lts[4], 40, ibc_e$e_40_lts[4], col=brew[4], lwd=1.5)
segments(40, ibc_e$e_40_lts[4], 45, ibc_e$e_45_lts[4], col=brew[4], lwd=1.5)
segments(45, ibc_e$e_45_lts[4], 50, ibc_e$e_50_lts[4], col=brew[4], lwd=1.5)
segments(50, ibc_e$e_50_lts[4], 55, ibc_e$e_55_lts[4], col=brew[4], lwd=1.5)
segments(55, ibc_e$e_55_lts[4], 60, ibc_e$e_60_lts[4], col=brew[4], lwd=1.5)
#PAG
lines(5, ibc_e$e_05_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(10, ibc_e$e_10_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(15, ibc_e$e_15_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(20, ibc_e$e_20_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(25, ibc_e$e_25_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(30, ibc_e$e_30_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(35, ibc_e$e_35_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(40, ibc_e$e_40_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(45, ibc_e$e_45_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(50, ibc_e$e_50_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(55, ibc_e$e_55_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(60, ibc_e$e_60_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
segments(0, 0, 5, ibc_e$e_05_lts[5], col=brew[5], lwd=1.5)
segments(5,  ibc_e$e_05_lts[5], 10, ibc_e$e_10_lts[5], col=brew[5], lwd=1.5)
segments(10, ibc_e$e_10_lts[5], 15, ibc_e$e_15_lts[5], col=brew[5], lwd=1.5)
segments(15, ibc_e$e_15_lts[5], 20, ibc_e$e_20_lts[5], col=brew[5], lwd=1.5)
segments(20, ibc_e$e_20_lts[5], 25, ibc_e$e_25_lts[5], col=brew[5], lwd=1.5)
segments(25, ibc_e$e_25_lts[5], 30, ibc_e$e_30_lts[5], col=brew[5], lwd=1.5)
segments(30, ibc_e$e_30_lts[5], 35, ibc_e$e_35_lts[5], col=brew[5], lwd=1.5)
segments(35, ibc_e$e_35_lts[5], 40, ibc_e$e_40_lts[5], col=brew[5], lwd=1.5)
segments(40, ibc_e$e_40_lts[5], 45, ibc_e$e_45_lts[5], col=brew[5], lwd=1.5)
segments(45, ibc_e$e_45_lts[5], 50, ibc_e$e_50_lts[5], col=brew[5], lwd=1.5)
segments(50, ibc_e$e_50_lts[5], 55, ibc_e$e_55_lts[5], col=brew[5], lwd=1.5)
segments(55, ibc_e$e_55_lts[5], 60, ibc_e$e_60_lts[5], col=brew[5], lwd=1.5)
#SCMPO
lines(5, ibc_e$e_05_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(10, ibc_e$e_10_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(15, ibc_e$e_15_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(20, ibc_e$e_20_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(25, ibc_e$e_25_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(30, ibc_e$e_30_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(35, ibc_e$e_35_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(40, ibc_e$e_40_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(45, ibc_e$e_45_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(50, ibc_e$e_50_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(55, ibc_e$e_55_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(60, ibc_e$e_60_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
segments(0, 0, 5, ibc_e$e_05_lts[6], col=brew[6], lwd=1.5)
segments(5,  ibc_e$e_05_lts[6], 10, ibc_e$e_10_lts[6], col=brew[6], lwd=1.5)
segments(10, ibc_e$e_10_lts[6], 15, ibc_e$e_15_lts[6], col=brew[6], lwd=1.5)
segments(15, ibc_e$e_15_lts[6], 20, ibc_e$e_20_lts[6], col=brew[6], lwd=1.5)
segments(20, ibc_e$e_20_lts[6], 25, ibc_e$e_25_lts[6], col=brew[6], lwd=1.5)
segments(25, ibc_e$e_25_lts[6], 30, ibc_e$e_30_lts[6], col=brew[6], lwd=1.5)
segments(30, ibc_e$e_30_lts[6], 35, ibc_e$e_35_lts[6], col=brew[6], lwd=1.5)
segments(35, ibc_e$e_35_lts[6], 40, ibc_e$e_40_lts[6], col=brew[6], lwd=1.5)
segments(40, ibc_e$e_40_lts[6], 45, ibc_e$e_45_lts[6], col=brew[6], lwd=1.5)
segments(45, ibc_e$e_45_lts[6], 50, ibc_e$e_50_lts[6], col=brew[6], lwd=1.5)
segments(50, ibc_e$e_50_lts[6], 55, ibc_e$e_55_lts[6], col=brew[6], lwd=1.5)
segments(55, ibc_e$e_55_lts[6], 60, ibc_e$e_60_lts[6], col=brew[6], lwd=1.5)
#SVMPO
lines(5, ibc_e$e_05_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(10, ibc_e$e_10_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(15, ibc_e$e_15_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(20, ibc_e$e_20_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(25, ibc_e$e_25_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(30, ibc_e$e_30_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(35, ibc_e$e_35_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(40, ibc_e$e_40_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(45, ibc_e$e_45_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(50, ibc_e$e_50_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(55, ibc_e$e_55_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(60, ibc_e$e_60_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
segments(0, 0, 5, ibc_e$e_05_lts[7], col=brew[7], lwd=1.5)
segments(5,  ibc_e$e_05_lts[7], 10, ibc_e$e_10_lts[7], col=brew[7], lwd=1.5)
segments(10, ibc_e$e_10_lts[7], 15, ibc_e$e_15_lts[7], col=brew[7], lwd=1.5)
segments(15, ibc_e$e_15_lts[7], 20, ibc_e$e_20_lts[7], col=brew[7], lwd=1.5)
segments(20, ibc_e$e_20_lts[7], 25, ibc_e$e_25_lts[7], col=brew[7], lwd=1.5)
segments(25, ibc_e$e_25_lts[7], 30, ibc_e$e_30_lts[7], col=brew[7], lwd=1.5)
segments(30, ibc_e$e_30_lts[7], 35, ibc_e$e_35_lts[7], col=brew[7], lwd=1.5)
segments(35, ibc_e$e_35_lts[7], 40, ibc_e$e_40_lts[7], col=brew[7], lwd=1.5)
segments(40, ibc_e$e_40_lts[7], 45, ibc_e$e_45_lts[7], col=brew[7], lwd=1.5)
segments(45, ibc_e$e_45_lts[7], 50, ibc_e$e_50_lts[7], col=brew[7], lwd=1.5)
segments(50, ibc_e$e_50_lts[7], 55, ibc_e$e_55_lts[7], col=brew[7], lwd=1.5)
segments(55, ibc_e$e_55_lts[7], 60, ibc_e$e_60_lts[7], col=brew[7], lwd=1.5)
#YMPO
lines(5, ibc_e$e_05_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(10, ibc_e$e_10_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(15, ibc_e$e_15_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(20, ibc_e$e_20_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(25, ibc_e$e_25_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(30, ibc_e$e_30_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(35, ibc_e$e_35_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(40, ibc_e$e_40_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(45, ibc_e$e_45_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(50, ibc_e$e_50_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(55, ibc_e$e_55_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(60, ibc_e$e_60_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
segments(0, 0, 5, ibc_e$e_05_lts[8], col=brew[8], lwd=1.5)
segments(5,  ibc_e$e_05_lts[8], 10, ibc_e$e_10_lts[8], col=brew[8], lwd=1.5)
segments(10, ibc_e$e_10_lts[8], 15, ibc_e$e_15_lts[8], col=brew[8], lwd=1.5)
segments(15, ibc_e$e_15_lts[8], 20, ibc_e$e_20_lts[8], col=brew[8], lwd=1.5)
segments(20, ibc_e$e_20_lts[8], 25, ibc_e$e_25_lts[8], col=brew[8], lwd=1.5)
segments(25, ibc_e$e_25_lts[8], 30, ibc_e$e_30_lts[8], col=brew[8], lwd=1.5)
segments(30, ibc_e$e_30_lts[8], 35, ibc_e$e_35_lts[8], col=brew[8], lwd=1.5)
segments(35, ibc_e$e_35_lts[8], 40, ibc_e$e_40_lts[8], col=brew[8], lwd=1.5)
segments(40, ibc_e$e_40_lts[8], 45, ibc_e$e_45_lts[8], col=brew[8], lwd=1.5)
segments(45, ibc_e$e_45_lts[8], 50, ibc_e$e_50_lts[8], col=brew[8], lwd=1.5)
segments(50, ibc_e$e_50_lts[8], 55, ibc_e$e_55_lts[8], col=brew[8], lwd=1.5)
segments(55, ibc_e$e_55_lts[8], 60, ibc_e$e_60_lts[8], col=brew[8], lwd=1.5)
# Create x-axis
axis(side=1, at=c(0,5,10,15,20,25,30,35,40,45,50,55,60), tick=TRUE, cex=1.2)
# Create legend
text(0, 50000, "Legend: IBC Cyclist", adj=0, cex=0.8, font=2, col="black")
points(0, 48000, col=brew[1], cex=0.8, pch=20)
segments(0, 48000, 5, 48000, col=brew[1], lty=1, lwd=1.5)
text(7, 48000, "CYMPO", adj=0, cex=0.8, col="black")
points(0, 47000, col=brew[3], cex=0.8, pch=20)
segments(0, 47000, 5, 47000, col=brew[3], lty=1, lwd=1.5)
text(7, 47000, "LHMPO", adj=0, cex=0.8, col="black")
points(0, 46000, col=brew[4], cex=0.8, pch=20)
segments(0, 46000, 5, 46000, col=brew[4], lty=1, lwd=1.5)
text(7, 46000, "MAG", adj=0, cex=0.8, col="black")
points(0, 45000, col=brew[2], cex=0.8, pch=20)
segments(0, 45000, 5, 45000, col=brew[2], lty=1, lwd=1.5)
text(7, 45000, "MetroPlan", adj=0, cex=0.8, col="black")
points(0, 44000, col=brew[5], cex=0.8, pch=20)
segments(0, 44000, 5, 44000, col=brew[5], lty=1, lwd=1.5)
text(7, 44000, "PAG", adj=0, cex=0.8, col="black")
points(0, 43000, col=brew[6], cex=0.8, pch=20)
segments(0, 43000, 5, 43000, col=brew[6], lty=1, lwd=1.5)
text(7, 43000, "SCMPO", adj=0, cex=0.8, col="black")
points(0, 42000, col=brew[7], cex=0.8, pch=20)
segments(0, 42000, 5, 42000, col=brew[7], lty=1, lwd=1.5)
text(7, 42000, "SVMPO", adj=0, cex=0.8, col="black")
points(0, 41000, col=brew[8], cex=0.8, pch=20)
segments(0, 41000, 5, 41000, col=brew[8], lty=1, lwd=1.5)
text(7, 41000, "YMPO", adj=0, cex=0.8, col="black")
rect(-1, 40000, 14, 51000, border="black", col=NA)
graphics.off()

### EAC Cyclist Type

# Graphic set-up
png(file=paste(dat_dir2, "plot_eac-employment.png", sep="/_analysis/_plots/"), bg="white", width=9, height=9, units='in', res=750)
plot(c(0,60), c(0,50000), type="n", xlab="Travel time (minutes)", ylab="Job accessibility", xaxt="n", cex=1.2)
# Plot points and lines
# CYMPO
lines(5, eac_e$e_05_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(10, eac_e$e_10_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(15, eac_e$e_15_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(20, eac_e$e_20_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(25, eac_e$e_25_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(30, eac_e$e_30_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(35, eac_e$e_35_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(40, eac_e$e_40_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(45, eac_e$e_45_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(50, eac_e$e_50_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(55, eac_e$e_55_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(60, eac_e$e_60_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
segments(0, 0, 5, eac_e$e_05_lts[1], col=brew[1], lwd=1.5)
segments(5,  eac_e$e_05_lts[1], 10, eac_e$e_10_lts[1], col=brew[1], lwd=1.5)
segments(10, eac_e$e_10_lts[1], 15, eac_e$e_15_lts[1], col=brew[1], lwd=1.5)
segments(15, eac_e$e_15_lts[1], 20, eac_e$e_20_lts[1], col=brew[1], lwd=1.5)
segments(20, eac_e$e_20_lts[1], 25, eac_e$e_25_lts[1], col=brew[1], lwd=1.5)
segments(25, eac_e$e_25_lts[1], 30, eac_e$e_30_lts[1], col=brew[1], lwd=1.5)
segments(30, eac_e$e_30_lts[1], 35, eac_e$e_35_lts[1], col=brew[1], lwd=1.5)
segments(35, eac_e$e_35_lts[1], 40, eac_e$e_40_lts[1], col=brew[1], lwd=1.5)
segments(40, eac_e$e_40_lts[1], 45, eac_e$e_45_lts[1], col=brew[1], lwd=1.5)
segments(45, eac_e$e_45_lts[1], 50, eac_e$e_50_lts[1], col=brew[1], lwd=1.5)
segments(50, eac_e$e_50_lts[1], 55, eac_e$e_55_lts[1], col=brew[1], lwd=1.5)
segments(55, eac_e$e_55_lts[1], 60, eac_e$e_60_lts[1], col=brew[1], lwd=1.5)
#FMPO
lines(5, eac_e$e_05_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(10, eac_e$e_10_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(15, eac_e$e_15_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(20, eac_e$e_20_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(25, eac_e$e_25_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(30, eac_e$e_30_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(35, eac_e$e_35_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(40, eac_e$e_40_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(45, eac_e$e_45_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(50, eac_e$e_50_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(55, eac_e$e_55_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(60, eac_e$e_60_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
segments(0, 0, 5, eac_e$e_05_lts[2], col=brew[2], lwd=1.5)
segments(5,  eac_e$e_05_lts[2], 10, eac_e$e_10_lts[2], col=brew[2], lwd=1.5)
segments(10, eac_e$e_10_lts[2], 15, eac_e$e_15_lts[2], col=brew[2], lwd=1.5)
segments(15, eac_e$e_15_lts[2], 20, eac_e$e_20_lts[2], col=brew[2], lwd=1.5)
segments(20, eac_e$e_20_lts[2], 25, eac_e$e_25_lts[2], col=brew[2], lwd=1.5)
segments(25, eac_e$e_25_lts[2], 30, eac_e$e_30_lts[2], col=brew[2], lwd=1.5)
segments(30, eac_e$e_30_lts[2], 35, eac_e$e_35_lts[2], col=brew[2], lwd=1.5)
segments(35, eac_e$e_35_lts[2], 40, eac_e$e_40_lts[2], col=brew[2], lwd=1.5)
segments(40, eac_e$e_40_lts[2], 45, eac_e$e_45_lts[2], col=brew[2], lwd=1.5)
segments(45, eac_e$e_45_lts[2], 50, eac_e$e_50_lts[2], col=brew[2], lwd=1.5)
segments(50, eac_e$e_50_lts[2], 55, eac_e$e_55_lts[2], col=brew[2], lwd=1.5)
segments(55, eac_e$e_55_lts[2], 60, eac_e$e_60_lts[2], col=brew[2], lwd=1.5)
#LHMPO
lines(5, eac_e$e_05_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(10, eac_e$e_10_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(15, eac_e$e_15_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(20, eac_e$e_20_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(25, eac_e$e_25_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(30, eac_e$e_30_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(35, eac_e$e_35_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(40, eac_e$e_40_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(45, eac_e$e_45_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(50, eac_e$e_50_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(55, eac_e$e_55_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(60, eac_e$e_60_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
segments(0, 0, 5, eac_e$e_05_lts[3], col=brew[3], lwd=1.5)
segments(5,  eac_e$e_05_lts[3], 10, eac_e$e_10_lts[3], col=brew[3], lwd=1.5)
segments(10, eac_e$e_10_lts[3], 15, eac_e$e_15_lts[3], col=brew[3], lwd=1.5)
segments(15, eac_e$e_15_lts[3], 20, eac_e$e_20_lts[3], col=brew[3], lwd=1.5)
segments(20, eac_e$e_20_lts[3], 25, eac_e$e_25_lts[3], col=brew[3], lwd=1.5)
segments(25, eac_e$e_25_lts[3], 30, eac_e$e_30_lts[3], col=brew[3], lwd=1.5)
segments(30, eac_e$e_30_lts[3], 35, eac_e$e_35_lts[3], col=brew[3], lwd=1.5)
segments(35, eac_e$e_35_lts[3], 40, eac_e$e_40_lts[3], col=brew[3], lwd=1.5)
segments(40, eac_e$e_40_lts[3], 45, eac_e$e_45_lts[3], col=brew[3], lwd=1.5)
segments(45, eac_e$e_45_lts[3], 50, eac_e$e_50_lts[3], col=brew[3], lwd=1.5)
segments(50, eac_e$e_50_lts[3], 55, eac_e$e_55_lts[3], col=brew[3], lwd=1.5)
segments(55, eac_e$e_55_lts[3], 60, eac_e$e_60_lts[3], col=brew[3], lwd=1.5)
#MAG
lines(5, eac_e$e_05_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(10, eac_e$e_10_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(15, eac_e$e_15_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(20, eac_e$e_20_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(25, eac_e$e_25_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(30, eac_e$e_30_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(35, eac_e$e_35_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(40, eac_e$e_40_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(45, eac_e$e_45_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(50, eac_e$e_50_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(55, eac_e$e_55_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(60, eac_e$e_60_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
segments(0, 0, 5, eac_e$e_05_lts[4], col=brew[4], lwd=1.5)
segments(5,  eac_e$e_05_lts[4], 10, eac_e$e_10_lts[4], col=brew[4], lwd=1.5)
segments(10, eac_e$e_10_lts[4], 15, eac_e$e_15_lts[4], col=brew[4], lwd=1.5)
segments(15, eac_e$e_15_lts[4], 20, eac_e$e_20_lts[4], col=brew[4], lwd=1.5)
segments(20, eac_e$e_20_lts[4], 25, eac_e$e_25_lts[4], col=brew[4], lwd=1.5)
segments(25, eac_e$e_25_lts[4], 30, eac_e$e_30_lts[4], col=brew[4], lwd=1.5)
segments(30, eac_e$e_30_lts[4], 35, eac_e$e_35_lts[4], col=brew[4], lwd=1.5)
segments(35, eac_e$e_35_lts[4], 40, eac_e$e_40_lts[4], col=brew[4], lwd=1.5)
segments(40, eac_e$e_40_lts[4], 45, eac_e$e_45_lts[4], col=brew[4], lwd=1.5)
segments(45, eac_e$e_45_lts[4], 50, eac_e$e_50_lts[4], col=brew[4], lwd=1.5)
segments(50, eac_e$e_50_lts[4], 55, eac_e$e_55_lts[4], col=brew[4], lwd=1.5)
segments(55, eac_e$e_55_lts[4], 60, eac_e$e_60_lts[4], col=brew[4], lwd=1.5)
#PAG
lines(5, eac_e$e_05_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(10, eac_e$e_10_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(15, eac_e$e_15_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(20, eac_e$e_20_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(25, eac_e$e_25_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(30, eac_e$e_30_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(35, eac_e$e_35_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(40, eac_e$e_40_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(45, eac_e$e_45_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(50, eac_e$e_50_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(55, eac_e$e_55_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(60, eac_e$e_60_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
segments(0, 0, 5, eac_e$e_05_lts[5], col=brew[5], lwd=1.5)
segments(5,  eac_e$e_05_lts[5], 10, eac_e$e_10_lts[5], col=brew[5], lwd=1.5)
segments(10, eac_e$e_10_lts[5], 15, eac_e$e_15_lts[5], col=brew[5], lwd=1.5)
segments(15, eac_e$e_15_lts[5], 20, eac_e$e_20_lts[5], col=brew[5], lwd=1.5)
segments(20, eac_e$e_20_lts[5], 25, eac_e$e_25_lts[5], col=brew[5], lwd=1.5)
segments(25, eac_e$e_25_lts[5], 30, eac_e$e_30_lts[5], col=brew[5], lwd=1.5)
segments(30, eac_e$e_30_lts[5], 35, eac_e$e_35_lts[5], col=brew[5], lwd=1.5)
segments(35, eac_e$e_35_lts[5], 40, eac_e$e_40_lts[5], col=brew[5], lwd=1.5)
segments(40, eac_e$e_40_lts[5], 45, eac_e$e_45_lts[5], col=brew[5], lwd=1.5)
segments(45, eac_e$e_45_lts[5], 50, eac_e$e_50_lts[5], col=brew[5], lwd=1.5)
segments(50, eac_e$e_50_lts[5], 55, eac_e$e_55_lts[5], col=brew[5], lwd=1.5)
segments(55, eac_e$e_55_lts[5], 60, eac_e$e_60_lts[5], col=brew[5], lwd=1.5)
#SCMPO
lines(5, eac_e$e_05_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(10, eac_e$e_10_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(15, eac_e$e_15_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(20, eac_e$e_20_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(25, eac_e$e_25_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(30, eac_e$e_30_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(35, eac_e$e_35_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(40, eac_e$e_40_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(45, eac_e$e_45_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(50, eac_e$e_50_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(55, eac_e$e_55_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(60, eac_e$e_60_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
segments(0, 0, 5, eac_e$e_05_lts[6], col=brew[6], lwd=1.5)
segments(5,  eac_e$e_05_lts[6], 10, eac_e$e_10_lts[6], col=brew[6], lwd=1.5)
segments(10, eac_e$e_10_lts[6], 15, eac_e$e_15_lts[6], col=brew[6], lwd=1.5)
segments(15, eac_e$e_15_lts[6], 20, eac_e$e_20_lts[6], col=brew[6], lwd=1.5)
segments(20, eac_e$e_20_lts[6], 25, eac_e$e_25_lts[6], col=brew[6], lwd=1.5)
segments(25, eac_e$e_25_lts[6], 30, eac_e$e_30_lts[6], col=brew[6], lwd=1.5)
segments(30, eac_e$e_30_lts[6], 35, eac_e$e_35_lts[6], col=brew[6], lwd=1.5)
segments(35, eac_e$e_35_lts[6], 40, eac_e$e_40_lts[6], col=brew[6], lwd=1.5)
segments(40, eac_e$e_40_lts[6], 45, eac_e$e_45_lts[6], col=brew[6], lwd=1.5)
segments(45, eac_e$e_45_lts[6], 50, eac_e$e_50_lts[6], col=brew[6], lwd=1.5)
segments(50, eac_e$e_50_lts[6], 55, eac_e$e_55_lts[6], col=brew[6], lwd=1.5)
segments(55, eac_e$e_55_lts[6], 60, eac_e$e_60_lts[6], col=brew[6], lwd=1.5)
#SVMPO
lines(5, eac_e$e_05_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(10, eac_e$e_10_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(15, eac_e$e_15_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(20, eac_e$e_20_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(25, eac_e$e_25_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(30, eac_e$e_30_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(35, eac_e$e_35_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(40, eac_e$e_40_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(45, eac_e$e_45_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(50, eac_e$e_50_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(55, eac_e$e_55_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(60, eac_e$e_60_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
segments(0, 0, 5, eac_e$e_05_lts[7], col=brew[7], lwd=1.5)
segments(5,  eac_e$e_05_lts[7], 10, eac_e$e_10_lts[7], col=brew[7], lwd=1.5)
segments(10, eac_e$e_10_lts[7], 15, eac_e$e_15_lts[7], col=brew[7], lwd=1.5)
segments(15, eac_e$e_15_lts[7], 20, eac_e$e_20_lts[7], col=brew[7], lwd=1.5)
segments(20, eac_e$e_20_lts[7], 25, eac_e$e_25_lts[7], col=brew[7], lwd=1.5)
segments(25, eac_e$e_25_lts[7], 30, eac_e$e_30_lts[7], col=brew[7], lwd=1.5)
segments(30, eac_e$e_30_lts[7], 35, eac_e$e_35_lts[7], col=brew[7], lwd=1.5)
segments(35, eac_e$e_35_lts[7], 40, eac_e$e_40_lts[7], col=brew[7], lwd=1.5)
segments(40, eac_e$e_40_lts[7], 45, eac_e$e_45_lts[7], col=brew[7], lwd=1.5)
segments(45, eac_e$e_45_lts[7], 50, eac_e$e_50_lts[7], col=brew[7], lwd=1.5)
segments(50, eac_e$e_50_lts[7], 55, eac_e$e_55_lts[7], col=brew[7], lwd=1.5)
segments(55, eac_e$e_55_lts[7], 60, eac_e$e_60_lts[7], col=brew[7], lwd=1.5)
#YMPO
lines(5, eac_e$e_05_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(10, eac_e$e_10_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(15, eac_e$e_15_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(20, eac_e$e_20_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(25, eac_e$e_25_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(30, eac_e$e_30_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(35, eac_e$e_35_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(40, eac_e$e_40_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(45, eac_e$e_45_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(50, eac_e$e_50_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(55, eac_e$e_55_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(60, eac_e$e_60_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
segments(0, 0, 5, eac_e$e_05_lts[8], col=brew[8], lwd=1.5)
segments(5,  eac_e$e_05_lts[8], 10, eac_e$e_10_lts[8], col=brew[8], lwd=1.5)
segments(10, eac_e$e_10_lts[8], 15, eac_e$e_15_lts[8], col=brew[8], lwd=1.5)
segments(15, eac_e$e_15_lts[8], 20, eac_e$e_20_lts[8], col=brew[8], lwd=1.5)
segments(20, eac_e$e_20_lts[8], 25, eac_e$e_25_lts[8], col=brew[8], lwd=1.5)
segments(25, eac_e$e_25_lts[8], 30, eac_e$e_30_lts[8], col=brew[8], lwd=1.5)
segments(30, eac_e$e_30_lts[8], 35, eac_e$e_35_lts[8], col=brew[8], lwd=1.5)
segments(35, eac_e$e_35_lts[8], 40, eac_e$e_40_lts[8], col=brew[8], lwd=1.5)
segments(40, eac_e$e_40_lts[8], 45, eac_e$e_45_lts[8], col=brew[8], lwd=1.5)
segments(45, eac_e$e_45_lts[8], 50, eac_e$e_50_lts[8], col=brew[8], lwd=1.5)
segments(50, eac_e$e_50_lts[8], 55, eac_e$e_55_lts[8], col=brew[8], lwd=1.5)
segments(55, eac_e$e_55_lts[8], 60, eac_e$e_60_lts[8], col=brew[8], lwd=1.5)
# Create x-axis
axis(side=1, at=c(0,5,10,15,20,25,30,35,40,45,50,55,60), tick=TRUE, cex=1.2)
# Create legend
text(0, 50000, "Legend: EAC Cyclist", adj=0, cex=0.8, font=2, col="black")
points(0, 48000, col=brew[1], cex=0.8, pch=20)
segments(0, 48000, 5, 48000, col=brew[1], lty=1, lwd=1.5)
text(7, 48000, "CYMPO", adj=0, cex=0.8, col="black")
points(0, 47000, col=brew[3], cex=0.8, pch=20)
segments(0, 47000, 5, 47000, col=brew[3], lty=1, lwd=1.5)
text(7, 47000, "LHMPO", adj=0, cex=0.8, col="black")
points(0, 46000, col=brew[4], cex=0.8, pch=20)
segments(0, 46000, 5, 46000, col=brew[4], lty=1, lwd=1.5)
text(7, 46000, "MAG", adj=0, cex=0.8, col="black")
points(0, 45000, col=brew[2], cex=0.8, pch=20)
segments(0, 45000, 5, 45000, col=brew[2], lty=1, lwd=1.5)
text(7, 45000, "MetroPlan", adj=0, cex=0.8, col="black")
points(0, 44000, col=brew[5], cex=0.8, pch=20)
segments(0, 44000, 5, 44000, col=brew[5], lty=1, lwd=1.5)
text(7, 44000, "PAG", adj=0, cex=0.8, col="black")
points(0, 43000, col=brew[6], cex=0.8, pch=20)
segments(0, 43000, 5, 43000, col=brew[6], lty=1, lwd=1.5)
text(7, 43000, "SCMPO", adj=0, cex=0.8, col="black")
points(0, 42000, col=brew[7], cex=0.8, pch=20)
segments(0, 42000, 5, 42000, col=brew[7], lty=1, lwd=1.5)
text(7, 42000, "SVMPO", adj=0, cex=0.8, col="black")
points(0, 41000, col=brew[8], cex=0.8, pch=20)
segments(0, 41000, 5, 41000, col=brew[8], lty=1, lwd=1.5)
text(7, 41000, "YMPO", adj=0, cex=0.8, col="black")
rect(-1, 40000, 14, 51000, border="black", col=NA)
graphics.off()

### SAF Cyclist Type

# Graphic set-up
png(file=paste(dat_dir2, "plot_saf-employment.png", sep="/_analysis/_plots/"), bg="white", width=9, height=9, units='in', res=750)
plot(c(0,60), c(0,50000), type="n", xlab="Travel time (minutes)", ylab="Job accessibility", xaxt="n", cex=1.2)
# Plot points and lines
# CYMPO
lines(5, saf_e$e_05_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(10, saf_e$e_10_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(15, saf_e$e_15_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(20, saf_e$e_20_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(25, saf_e$e_25_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(30, saf_e$e_30_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(35, saf_e$e_35_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(40, saf_e$e_40_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(45, saf_e$e_45_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(50, saf_e$e_50_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(55, saf_e$e_55_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(60, saf_e$e_60_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
segments(0, 0, 5, saf_e$e_05_lts[1], col=brew[1], lwd=1.5)
segments(5,  saf_e$e_05_lts[1], 10, saf_e$e_10_lts[1], col=brew[1], lwd=1.5)
segments(10, saf_e$e_10_lts[1], 15, saf_e$e_15_lts[1], col=brew[1], lwd=1.5)
segments(15, saf_e$e_15_lts[1], 20, saf_e$e_20_lts[1], col=brew[1], lwd=1.5)
segments(20, saf_e$e_20_lts[1], 25, saf_e$e_25_lts[1], col=brew[1], lwd=1.5)
segments(25, saf_e$e_25_lts[1], 30, saf_e$e_30_lts[1], col=brew[1], lwd=1.5)
segments(30, saf_e$e_30_lts[1], 35, saf_e$e_35_lts[1], col=brew[1], lwd=1.5)
segments(35, saf_e$e_35_lts[1], 40, saf_e$e_40_lts[1], col=brew[1], lwd=1.5)
segments(40, saf_e$e_40_lts[1], 45, saf_e$e_45_lts[1], col=brew[1], lwd=1.5)
segments(45, saf_e$e_45_lts[1], 50, saf_e$e_50_lts[1], col=brew[1], lwd=1.5)
segments(50, saf_e$e_50_lts[1], 55, saf_e$e_55_lts[1], col=brew[1], lwd=1.5)
segments(55, saf_e$e_55_lts[1], 60, saf_e$e_60_lts[1], col=brew[1], lwd=1.5)
#FMPO
lines(5, saf_e$e_05_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(10, saf_e$e_10_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(15, saf_e$e_15_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(20, saf_e$e_20_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(25, saf_e$e_25_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(30, saf_e$e_30_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(35, saf_e$e_35_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(40, saf_e$e_40_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(45, saf_e$e_45_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(50, saf_e$e_50_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(55, saf_e$e_55_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(60, saf_e$e_60_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
segments(0, 0, 5, saf_e$e_05_lts[2], col=brew[2], lwd=1.5)
segments(5,  saf_e$e_05_lts[2], 10, saf_e$e_10_lts[2], col=brew[2], lwd=1.5)
segments(10, saf_e$e_10_lts[2], 15, saf_e$e_15_lts[2], col=brew[2], lwd=1.5)
segments(15, saf_e$e_15_lts[2], 20, saf_e$e_20_lts[2], col=brew[2], lwd=1.5)
segments(20, saf_e$e_20_lts[2], 25, saf_e$e_25_lts[2], col=brew[2], lwd=1.5)
segments(25, saf_e$e_25_lts[2], 30, saf_e$e_30_lts[2], col=brew[2], lwd=1.5)
segments(30, saf_e$e_30_lts[2], 35, saf_e$e_35_lts[2], col=brew[2], lwd=1.5)
segments(35, saf_e$e_35_lts[2], 40, saf_e$e_40_lts[2], col=brew[2], lwd=1.5)
segments(40, saf_e$e_40_lts[2], 45, saf_e$e_45_lts[2], col=brew[2], lwd=1.5)
segments(45, saf_e$e_45_lts[2], 50, saf_e$e_50_lts[2], col=brew[2], lwd=1.5)
segments(50, saf_e$e_50_lts[2], 55, saf_e$e_55_lts[2], col=brew[2], lwd=1.5)
segments(55, saf_e$e_55_lts[2], 60, saf_e$e_60_lts[2], col=brew[2], lwd=1.5)
#LHMPO
lines(5, saf_e$e_05_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(10, saf_e$e_10_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(15, saf_e$e_15_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(20, saf_e$e_20_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(25, saf_e$e_25_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(30, saf_e$e_30_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(35, saf_e$e_35_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(40, saf_e$e_40_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(45, saf_e$e_45_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(50, saf_e$e_50_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(55, saf_e$e_55_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(60, saf_e$e_60_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
segments(0, 0, 5, saf_e$e_05_lts[3], col=brew[3], lwd=1.5)
segments(5,  saf_e$e_05_lts[3], 10, saf_e$e_10_lts[3], col=brew[3], lwd=1.5)
segments(10, saf_e$e_10_lts[3], 15, saf_e$e_15_lts[3], col=brew[3], lwd=1.5)
segments(15, saf_e$e_15_lts[3], 20, saf_e$e_20_lts[3], col=brew[3], lwd=1.5)
segments(20, saf_e$e_20_lts[3], 25, saf_e$e_25_lts[3], col=brew[3], lwd=1.5)
segments(25, saf_e$e_25_lts[3], 30, saf_e$e_30_lts[3], col=brew[3], lwd=1.5)
segments(30, saf_e$e_30_lts[3], 35, saf_e$e_35_lts[3], col=brew[3], lwd=1.5)
segments(35, saf_e$e_35_lts[3], 40, saf_e$e_40_lts[3], col=brew[3], lwd=1.5)
segments(40, saf_e$e_40_lts[3], 45, saf_e$e_45_lts[3], col=brew[3], lwd=1.5)
segments(45, saf_e$e_45_lts[3], 50, saf_e$e_50_lts[3], col=brew[3], lwd=1.5)
segments(50, saf_e$e_50_lts[3], 55, saf_e$e_55_lts[3], col=brew[3], lwd=1.5)
segments(55, saf_e$e_55_lts[3], 60, saf_e$e_60_lts[3], col=brew[3], lwd=1.5)
#MAG
lines(5, saf_e$e_05_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(10, saf_e$e_10_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(15, saf_e$e_15_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(20, saf_e$e_20_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(25, saf_e$e_25_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(30, saf_e$e_30_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(35, saf_e$e_35_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(40, saf_e$e_40_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(45, saf_e$e_45_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(50, saf_e$e_50_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(55, saf_e$e_55_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(60, saf_e$e_60_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
segments(0, 0, 5, saf_e$e_05_lts[4], col=brew[4], lwd=1.5)
segments(5,  saf_e$e_05_lts[4], 10, saf_e$e_10_lts[4], col=brew[4], lwd=1.5)
segments(10, saf_e$e_10_lts[4], 15, saf_e$e_15_lts[4], col=brew[4], lwd=1.5)
segments(15, saf_e$e_15_lts[4], 20, saf_e$e_20_lts[4], col=brew[4], lwd=1.5)
segments(20, saf_e$e_20_lts[4], 25, saf_e$e_25_lts[4], col=brew[4], lwd=1.5)
segments(25, saf_e$e_25_lts[4], 30, saf_e$e_30_lts[4], col=brew[4], lwd=1.5)
segments(30, saf_e$e_30_lts[4], 35, saf_e$e_35_lts[4], col=brew[4], lwd=1.5)
segments(35, saf_e$e_35_lts[4], 40, saf_e$e_40_lts[4], col=brew[4], lwd=1.5)
segments(40, saf_e$e_40_lts[4], 45, saf_e$e_45_lts[4], col=brew[4], lwd=1.5)
segments(45, saf_e$e_45_lts[4], 50, saf_e$e_50_lts[4], col=brew[4], lwd=1.5)
segments(50, saf_e$e_50_lts[4], 55, saf_e$e_55_lts[4], col=brew[4], lwd=1.5)
segments(55, saf_e$e_55_lts[4], 60, saf_e$e_60_lts[4], col=brew[4], lwd=1.5)
#PAG
lines(5, saf_e$e_05_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(10, saf_e$e_10_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(15, saf_e$e_15_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(20, saf_e$e_20_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(25, saf_e$e_25_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(30, saf_e$e_30_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(35, saf_e$e_35_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(40, saf_e$e_40_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(45, saf_e$e_45_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(50, saf_e$e_50_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(55, saf_e$e_55_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(60, saf_e$e_60_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
segments(0, 0, 5, saf_e$e_05_lts[5], col=brew[5], lwd=1.5)
segments(5,  saf_e$e_05_lts[5], 10, saf_e$e_10_lts[5], col=brew[5], lwd=1.5)
segments(10, saf_e$e_10_lts[5], 15, saf_e$e_15_lts[5], col=brew[5], lwd=1.5)
segments(15, saf_e$e_15_lts[5], 20, saf_e$e_20_lts[5], col=brew[5], lwd=1.5)
segments(20, saf_e$e_20_lts[5], 25, saf_e$e_25_lts[5], col=brew[5], lwd=1.5)
segments(25, saf_e$e_25_lts[5], 30, saf_e$e_30_lts[5], col=brew[5], lwd=1.5)
segments(30, saf_e$e_30_lts[5], 35, saf_e$e_35_lts[5], col=brew[5], lwd=1.5)
segments(35, saf_e$e_35_lts[5], 40, saf_e$e_40_lts[5], col=brew[5], lwd=1.5)
segments(40, saf_e$e_40_lts[5], 45, saf_e$e_45_lts[5], col=brew[5], lwd=1.5)
segments(45, saf_e$e_45_lts[5], 50, saf_e$e_50_lts[5], col=brew[5], lwd=1.5)
segments(50, saf_e$e_50_lts[5], 55, saf_e$e_55_lts[5], col=brew[5], lwd=1.5)
segments(55, saf_e$e_55_lts[5], 60, saf_e$e_60_lts[5], col=brew[5], lwd=1.5)
#SCMPO
lines(5, saf_e$e_05_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(10, saf_e$e_10_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(15, saf_e$e_15_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(20, saf_e$e_20_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(25, saf_e$e_25_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(30, saf_e$e_30_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(35, saf_e$e_35_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(40, saf_e$e_40_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(45, saf_e$e_45_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(50, saf_e$e_50_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(55, saf_e$e_55_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(60, saf_e$e_60_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
segments(0, 0, 5, saf_e$e_05_lts[6], col=brew[6], lwd=1.5)
segments(5,  saf_e$e_05_lts[6], 10, saf_e$e_10_lts[6], col=brew[6], lwd=1.5)
segments(10, saf_e$e_10_lts[6], 15, saf_e$e_15_lts[6], col=brew[6], lwd=1.5)
segments(15, saf_e$e_15_lts[6], 20, saf_e$e_20_lts[6], col=brew[6], lwd=1.5)
segments(20, saf_e$e_20_lts[6], 25, saf_e$e_25_lts[6], col=brew[6], lwd=1.5)
segments(25, saf_e$e_25_lts[6], 30, saf_e$e_30_lts[6], col=brew[6], lwd=1.5)
segments(30, saf_e$e_30_lts[6], 35, saf_e$e_35_lts[6], col=brew[6], lwd=1.5)
segments(35, saf_e$e_35_lts[6], 40, saf_e$e_40_lts[6], col=brew[6], lwd=1.5)
segments(40, saf_e$e_40_lts[6], 45, saf_e$e_45_lts[6], col=brew[6], lwd=1.5)
segments(45, saf_e$e_45_lts[6], 50, saf_e$e_50_lts[6], col=brew[6], lwd=1.5)
segments(50, saf_e$e_50_lts[6], 55, saf_e$e_55_lts[6], col=brew[6], lwd=1.5)
segments(55, saf_e$e_55_lts[6], 60, saf_e$e_60_lts[6], col=brew[6], lwd=1.5)
#SVMPO
lines(5, saf_e$e_05_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(10, saf_e$e_10_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(15, saf_e$e_15_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(20, saf_e$e_20_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(25, saf_e$e_25_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(30, saf_e$e_30_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(35, saf_e$e_35_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(40, saf_e$e_40_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(45, saf_e$e_45_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(50, saf_e$e_50_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(55, saf_e$e_55_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(60, saf_e$e_60_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
segments(0, 0, 5, saf_e$e_05_lts[7], col=brew[7], lwd=1.5)
segments(5,  saf_e$e_05_lts[7], 10, saf_e$e_10_lts[7], col=brew[7], lwd=1.5)
segments(10, saf_e$e_10_lts[7], 15, saf_e$e_15_lts[7], col=brew[7], lwd=1.5)
segments(15, saf_e$e_15_lts[7], 20, saf_e$e_20_lts[7], col=brew[7], lwd=1.5)
segments(20, saf_e$e_20_lts[7], 25, saf_e$e_25_lts[7], col=brew[7], lwd=1.5)
segments(25, saf_e$e_25_lts[7], 30, saf_e$e_30_lts[7], col=brew[7], lwd=1.5)
segments(30, saf_e$e_30_lts[7], 35, saf_e$e_35_lts[7], col=brew[7], lwd=1.5)
segments(35, saf_e$e_35_lts[7], 40, saf_e$e_40_lts[7], col=brew[7], lwd=1.5)
segments(40, saf_e$e_40_lts[7], 45, saf_e$e_45_lts[7], col=brew[7], lwd=1.5)
segments(45, saf_e$e_45_lts[7], 50, saf_e$e_50_lts[7], col=brew[7], lwd=1.5)
segments(50, saf_e$e_50_lts[7], 55, saf_e$e_55_lts[7], col=brew[7], lwd=1.5)
segments(55, saf_e$e_55_lts[7], 60, saf_e$e_60_lts[7], col=brew[7], lwd=1.5)
#YMPO
lines(5, saf_e$e_05_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(10, saf_e$e_10_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(15, saf_e$e_15_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(20, saf_e$e_20_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(25, saf_e$e_25_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(30, saf_e$e_30_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(35, saf_e$e_35_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(40, saf_e$e_40_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(45, saf_e$e_45_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(50, saf_e$e_50_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(55, saf_e$e_55_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(60, saf_e$e_60_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
segments(0, 0, 5, saf_e$e_05_lts[8], col=brew[8], lwd=1.5)
segments(5,  saf_e$e_05_lts[8], 10, saf_e$e_10_lts[8], col=brew[8], lwd=1.5)
segments(10, saf_e$e_10_lts[8], 15, saf_e$e_15_lts[8], col=brew[8], lwd=1.5)
segments(15, saf_e$e_15_lts[8], 20, saf_e$e_20_lts[8], col=brew[8], lwd=1.5)
segments(20, saf_e$e_20_lts[8], 25, saf_e$e_25_lts[8], col=brew[8], lwd=1.5)
segments(25, saf_e$e_25_lts[8], 30, saf_e$e_30_lts[8], col=brew[8], lwd=1.5)
segments(30, saf_e$e_30_lts[8], 35, saf_e$e_35_lts[8], col=brew[8], lwd=1.5)
segments(35, saf_e$e_35_lts[8], 40, saf_e$e_40_lts[8], col=brew[8], lwd=1.5)
segments(40, saf_e$e_40_lts[8], 45, saf_e$e_45_lts[8], col=brew[8], lwd=1.5)
segments(45, saf_e$e_45_lts[8], 50, saf_e$e_50_lts[8], col=brew[8], lwd=1.5)
segments(50, saf_e$e_50_lts[8], 55, saf_e$e_55_lts[8], col=brew[8], lwd=1.5)
segments(55, saf_e$e_55_lts[8], 60, saf_e$e_60_lts[8], col=brew[8], lwd=1.5)
# Create x-axis
axis(side=1, at=c(0,5,10,15,20,25,30,35,40,45,50,55,60), tick=TRUE, cex=1.2)
# Create legend
text(0, 50000, "Legend: SAF Cyclist", adj=0, cex=0.8, font=2, col="black")
points(0, 48000, col=brew[1], cex=0.8, pch=20)
segments(0, 48000, 5, 48000, col=brew[1], lty=1, lwd=1.5)
text(7, 48000, "CYMPO", adj=0, cex=0.8, col="black")
points(0, 47000, col=brew[3], cex=0.8, pch=20)
segments(0, 47000, 5, 47000, col=brew[3], lty=1, lwd=1.5)
text(7, 47000, "LHMPO", adj=0, cex=0.8, col="black")
points(0, 46000, col=brew[4], cex=0.8, pch=20)
segments(0, 46000, 5, 46000, col=brew[4], lty=1, lwd=1.5)
text(7, 46000, "MAG", adj=0, cex=0.8, col="black")
points(0, 45000, col=brew[2], cex=0.8, pch=20)
segments(0, 45000, 5, 45000, col=brew[2], lty=1, lwd=1.5)
text(7, 45000, "MetroPlan", adj=0, cex=0.8, col="black")
points(0, 44000, col=brew[5], cex=0.8, pch=20)
segments(0, 44000, 5, 44000, col=brew[5], lty=1, lwd=1.5)
text(7, 44000, "PAG", adj=0, cex=0.8, col="black")
points(0, 43000, col=brew[6], cex=0.8, pch=20)
segments(0, 43000, 5, 43000, col=brew[6], lty=1, lwd=1.5)
text(7, 43000, "SCMPO", adj=0, cex=0.8, col="black")
points(0, 42000, col=brew[7], cex=0.8, pch=20)
segments(0, 42000, 5, 42000, col=brew[7], lty=1, lwd=1.5)
text(7, 42000, "SVMPO", adj=0, cex=0.8, col="black")
points(0, 41000, col=brew[8], cex=0.8, pch=20)
segments(0, 41000, 5, 41000, col=brew[8], lty=1, lwd=1.5)
text(7, 41000, "YMPO", adj=0, cex=0.8, col="black")
rect(-1, 40000, 14, 51000, border="black", col=NA)
graphics.off()

####################
### MARKETS
####################

# Summarize grocery market access by cyclist type
ibc_m <- dat_ibc_m %>% group_by(mpo) %>%
  summarise(m_05=mean(m_00_05,na.rm=T),m_10=mean(m_05_10,na.rm=T),m_15=mean(m_10_15,na.rm=T),m_20=mean(m_15_20,na.rm=T),
            m_25=mean(m_20_25,na.rm=T),m_30=mean(m_25_30,na.rm=T),m_35=mean(m_30_35,na.rm=T),m_40=mean(m_35_40,na.rm=T),
            m_45=mean(m_40_45,na.rm=T),m_50=mean(m_45_50,na.rm=T),m_55=mean(m_50_55,na.rm=T),m_60=mean(m_55_60,na.rm=T),
            m_05_lts=mean(m_00_05_lt,na.rm=T),m_10_lts=mean(m_05_10_lt,na.rm=T),m_15_lts=mean(m_10_15_lt,na.rm=T),m_20_lts=mean(m_15_20_lt,na.rm=T),
            m_25_lts=mean(m_20_25_lt,na.rm=T),m_30_lts=mean(m_25_30_lt,na.rm=T),m_35_lts=mean(m_30_35_lt,na.rm=T),m_40_lts=mean(m_35_40_lt,na.rm=T),
            m_45_lts=mean(m_40_45_lt,na.rm=T),m_50_lts=mean(m_45_50_lt,na.rm=T),m_55_lts=mean(m_50_55_lt,na.rm=T),m_60_lts=mean(m_55_60_lt,na.rm=T))
eac_m <- dat_eac_m %>% group_by(mpo) %>%
  summarise(m_05=mean(m_00_05,na.rm=T),m_10=mean(m_05_10,na.rm=T),m_15=mean(m_10_15,na.rm=T),m_20=mean(m_15_20,na.rm=T),
            m_25=mean(m_20_25,na.rm=T),m_30=mean(m_25_30,na.rm=T),m_35=mean(m_30_35,na.rm=T),m_40=mean(m_35_40,na.rm=T),
            m_45=mean(m_40_45,na.rm=T),m_50=mean(m_45_50,na.rm=T),m_55=mean(m_50_55,na.rm=T),m_60=mean(m_55_60,na.rm=T),
            m_05_lts=mean(m_00_05_lt,na.rm=T),m_10_lts=mean(m_05_10_lt,na.rm=T),m_15_lts=mean(m_10_15_lt,na.rm=T),m_20_lts=mean(m_15_20_lt,na.rm=T),
            m_25_lts=mean(m_20_25_lt,na.rm=T),m_30_lts=mean(m_25_30_lt,na.rm=T),m_35_lts=mean(m_30_35_lt,na.rm=T),m_40_lts=mean(m_35_40_lt,na.rm=T),
            m_45_lts=mean(m_40_45_lt,na.rm=T),m_50_lts=mean(m_45_50_lt,na.rm=T),m_55_lts=mean(m_50_55_lt,na.rm=T),m_60_lts=mean(m_55_60_lt,na.rm=T))
saf_m <- dat_saf_m %>% group_by(mpo) %>%
  summarise(m_05=mean(m_00_05,na.rm=T),m_10=mean(m_05_10,na.rm=T),m_15=mean(m_10_15,na.rm=T),m_20=mean(m_15_20,na.rm=T),
            m_25=mean(m_20_25,na.rm=T),m_30=mean(m_25_30,na.rm=T),m_35=mean(m_30_35,na.rm=T),m_40=mean(m_35_40,na.rm=T),
            m_45=mean(m_40_45,na.rm=T),m_50=mean(m_45_50,na.rm=T),m_55=mean(m_50_55,na.rm=T),m_60=mean(m_55_60,na.rm=T),
            m_05_lts=mean(m_00_05_lt,na.rm=T),m_10_lts=mean(m_05_10_lt,na.rm=T),m_15_lts=mean(m_10_15_lt,na.rm=T),m_20_lts=mean(m_15_20_lt,na.rm=T),
            m_25_lts=mean(m_20_25_lt,na.rm=T),m_30_lts=mean(m_25_30_lt,na.rm=T),m_35_lts=mean(m_30_35_lt,na.rm=T),m_40_lts=mean(m_35_40_lt,na.rm=T),
            m_45_lts=mean(m_40_45_lt,na.rm=T),m_50_lts=mean(m_45_50_lt,na.rm=T),m_55_lts=mean(m_50_55_lt,na.rm=T),m_60_lts=mean(m_55_60_lt,na.rm=T))

### IBC Cyclist Type

# Graphic set-up
png(file=paste(dat_dir2, "plot_ibc-markets.png", sep="/_analysis/_plots/"), bg="white", width=9, height=9, units='in', res=750)
plot(c(0,60), c(0,10), type="n", xlab="Travel time (minutes)", ylab="Market accessibility", xaxt="n", cex=1.2)
# Plot points and lines
# CYMPO
lines(5, ibc_m$m_05_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(10, ibc_m$m_10_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(15, ibc_m$m_15_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(20, ibc_m$m_20_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(25, ibc_m$m_25_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(30, ibc_m$m_30_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(35, ibc_m$m_35_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(40, ibc_m$m_40_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(45, ibc_m$m_45_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(50, ibc_m$m_50_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(55, ibc_m$m_55_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(60, ibc_m$m_60_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
segments(0, 0, 5, ibc_m$m_05_lts[1], col=brew[1], lwd=1.5)
segments(5,  ibc_m$m_05_lts[1], 10, ibc_m$m_10_lts[1], col=brew[1], lwd=1.5)
segments(10, ibc_m$m_10_lts[1], 15, ibc_m$m_15_lts[1], col=brew[1], lwd=1.5)
segments(15, ibc_m$m_15_lts[1], 20, ibc_m$m_20_lts[1], col=brew[1], lwd=1.5)
segments(20, ibc_m$m_20_lts[1], 25, ibc_m$m_25_lts[1], col=brew[1], lwd=1.5)
segments(25, ibc_m$m_25_lts[1], 30, ibc_m$m_30_lts[1], col=brew[1], lwd=1.5)
segments(30, ibc_m$m_30_lts[1], 35, ibc_m$m_35_lts[1], col=brew[1], lwd=1.5)
segments(35, ibc_m$m_35_lts[1], 40, ibc_m$m_40_lts[1], col=brew[1], lwd=1.5)
segments(40, ibc_m$m_40_lts[1], 45, ibc_m$m_45_lts[1], col=brew[1], lwd=1.5)
segments(45, ibc_m$m_45_lts[1], 50, ibc_m$m_50_lts[1], col=brew[1], lwd=1.5)
segments(50, ibc_m$m_50_lts[1], 55, ibc_m$m_55_lts[1], col=brew[1], lwd=1.5)
segments(55, ibc_m$m_55_lts[1], 60, ibc_m$m_60_lts[1], col=brew[1], lwd=1.5)
#FMPO
lines(5, ibc_m$m_05_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(10, ibc_m$m_10_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(15, ibc_m$m_15_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(20, ibc_m$m_20_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(25, ibc_m$m_25_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(30, ibc_m$m_30_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(35, ibc_m$m_35_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(40, ibc_m$m_40_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(45, ibc_m$m_45_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(50, ibc_m$m_50_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(55, ibc_m$m_55_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(60, ibc_m$m_60_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
segments(0, 0, 5, ibc_m$m_05_lts[2], col=brew[2], lwd=1.5)
segments(5,  ibc_m$m_05_lts[2], 10, ibc_m$m_10_lts[2], col=brew[2], lwd=1.5)
segments(10, ibc_m$m_10_lts[2], 15, ibc_m$m_15_lts[2], col=brew[2], lwd=1.5)
segments(15, ibc_m$m_15_lts[2], 20, ibc_m$m_20_lts[2], col=brew[2], lwd=1.5)
segments(20, ibc_m$m_20_lts[2], 25, ibc_m$m_25_lts[2], col=brew[2], lwd=1.5)
segments(25, ibc_m$m_25_lts[2], 30, ibc_m$m_30_lts[2], col=brew[2], lwd=1.5)
segments(30, ibc_m$m_30_lts[2], 35, ibc_m$m_35_lts[2], col=brew[2], lwd=1.5)
segments(35, ibc_m$m_35_lts[2], 40, ibc_m$m_40_lts[2], col=brew[2], lwd=1.5)
segments(40, ibc_m$m_40_lts[2], 45, ibc_m$m_45_lts[2], col=brew[2], lwd=1.5)
segments(45, ibc_m$m_45_lts[2], 50, ibc_m$m_50_lts[2], col=brew[2], lwd=1.5)
segments(50, ibc_m$m_50_lts[2], 55, ibc_m$m_55_lts[2], col=brew[2], lwd=1.5)
segments(55, ibc_m$m_55_lts[2], 60, ibc_m$m_60_lts[2], col=brew[2], lwd=1.5)
#LHMPO
lines(5, ibc_m$m_05_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(10, ibc_m$m_10_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(15, ibc_m$m_15_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(20, ibc_m$m_20_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(25, ibc_m$m_25_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(30, ibc_m$m_30_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(35, ibc_m$m_35_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(40, ibc_m$m_40_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(45, ibc_m$m_45_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(50, ibc_m$m_50_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(55, ibc_m$m_55_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(60, ibc_m$m_60_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
segments(0, 0, 5, ibc_m$m_05_lts[3], col=brew[3], lwd=1.5)
segments(5,  ibc_m$m_05_lts[3], 10, ibc_m$m_10_lts[3], col=brew[3], lwd=1.5)
segments(10, ibc_m$m_10_lts[3], 15, ibc_m$m_15_lts[3], col=brew[3], lwd=1.5)
segments(15, ibc_m$m_15_lts[3], 20, ibc_m$m_20_lts[3], col=brew[3], lwd=1.5)
segments(20, ibc_m$m_20_lts[3], 25, ibc_m$m_25_lts[3], col=brew[3], lwd=1.5)
segments(25, ibc_m$m_25_lts[3], 30, ibc_m$m_30_lts[3], col=brew[3], lwd=1.5)
segments(30, ibc_m$m_30_lts[3], 35, ibc_m$m_35_lts[3], col=brew[3], lwd=1.5)
segments(35, ibc_m$m_35_lts[3], 40, ibc_m$m_40_lts[3], col=brew[3], lwd=1.5)
segments(40, ibc_m$m_40_lts[3], 45, ibc_m$m_45_lts[3], col=brew[3], lwd=1.5)
segments(45, ibc_m$m_45_lts[3], 50, ibc_m$m_50_lts[3], col=brew[3], lwd=1.5)
segments(50, ibc_m$m_50_lts[3], 55, ibc_m$m_55_lts[3], col=brew[3], lwd=1.5)
segments(55, ibc_m$m_55_lts[3], 60, ibc_m$m_60_lts[3], col=brew[3], lwd=1.5)
#MAG
lines(5, ibc_m$m_05_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(10, ibc_m$m_10_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(15, ibc_m$m_15_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(20, ibc_m$m_20_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(25, ibc_m$m_25_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(30, ibc_m$m_30_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(35, ibc_m$m_35_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(40, ibc_m$m_40_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(45, ibc_m$m_45_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(50, ibc_m$m_50_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(55, ibc_m$m_55_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(60, ibc_m$m_60_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
segments(0, 0, 5, ibc_m$m_05_lts[4], col=brew[4], lwd=1.5)
segments(5,  ibc_m$m_05_lts[4], 10, ibc_m$m_10_lts[4], col=brew[4], lwd=1.5)
segments(10, ibc_m$m_10_lts[4], 15, ibc_m$m_15_lts[4], col=brew[4], lwd=1.5)
segments(15, ibc_m$m_15_lts[4], 20, ibc_m$m_20_lts[4], col=brew[4], lwd=1.5)
segments(20, ibc_m$m_20_lts[4], 25, ibc_m$m_25_lts[4], col=brew[4], lwd=1.5)
segments(25, ibc_m$m_25_lts[4], 30, ibc_m$m_30_lts[4], col=brew[4], lwd=1.5)
segments(30, ibc_m$m_30_lts[4], 35, ibc_m$m_35_lts[4], col=brew[4], lwd=1.5)
segments(35, ibc_m$m_35_lts[4], 40, ibc_m$m_40_lts[4], col=brew[4], lwd=1.5)
segments(40, ibc_m$m_40_lts[4], 45, ibc_m$m_45_lts[4], col=brew[4], lwd=1.5)
segments(45, ibc_m$m_45_lts[4], 50, ibc_m$m_50_lts[4], col=brew[4], lwd=1.5)
segments(50, ibc_m$m_50_lts[4], 55, ibc_m$m_55_lts[4], col=brew[4], lwd=1.5)
segments(55, ibc_m$m_55_lts[4], 60, ibc_m$m_60_lts[4], col=brew[4], lwd=1.5)
#PAG
lines(5, ibc_m$m_05_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(10, ibc_m$m_10_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(15, ibc_m$m_15_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(20, ibc_m$m_20_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(25, ibc_m$m_25_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(30, ibc_m$m_30_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(35, ibc_m$m_35_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(40, ibc_m$m_40_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(45, ibc_m$m_45_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(50, ibc_m$m_50_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(55, ibc_m$m_55_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(60, ibc_m$m_60_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
segments(0, 0, 5, ibc_m$m_05_lts[5], col=brew[5], lwd=1.5)
segments(5,  ibc_m$m_05_lts[5], 10, ibc_m$m_10_lts[5], col=brew[5], lwd=1.5)
segments(10, ibc_m$m_10_lts[5], 15, ibc_m$m_15_lts[5], col=brew[5], lwd=1.5)
segments(15, ibc_m$m_15_lts[5], 20, ibc_m$m_20_lts[5], col=brew[5], lwd=1.5)
segments(20, ibc_m$m_20_lts[5], 25, ibc_m$m_25_lts[5], col=brew[5], lwd=1.5)
segments(25, ibc_m$m_25_lts[5], 30, ibc_m$m_30_lts[5], col=brew[5], lwd=1.5)
segments(30, ibc_m$m_30_lts[5], 35, ibc_m$m_35_lts[5], col=brew[5], lwd=1.5)
segments(35, ibc_m$m_35_lts[5], 40, ibc_m$m_40_lts[5], col=brew[5], lwd=1.5)
segments(40, ibc_m$m_40_lts[5], 45, ibc_m$m_45_lts[5], col=brew[5], lwd=1.5)
segments(45, ibc_m$m_45_lts[5], 50, ibc_m$m_50_lts[5], col=brew[5], lwd=1.5)
segments(50, ibc_m$m_50_lts[5], 55, ibc_m$m_55_lts[5], col=brew[5], lwd=1.5)
segments(55, ibc_m$m_55_lts[5], 60, ibc_m$m_60_lts[5], col=brew[5], lwd=1.5)
#SCMPO
lines(5, ibc_m$m_05_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(10, ibc_m$m_10_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(15, ibc_m$m_15_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(20, ibc_m$m_20_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(25, ibc_m$m_25_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(30, ibc_m$m_30_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(35, ibc_m$m_35_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(40, ibc_m$m_40_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(45, ibc_m$m_45_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(50, ibc_m$m_50_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(55, ibc_m$m_55_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(60, ibc_m$m_60_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
segments(0, 0, 5, ibc_m$m_05_lts[6], col=brew[6], lwd=1.5)
segments(5,  ibc_m$m_05_lts[6], 10, ibc_m$m_10_lts[6], col=brew[6], lwd=1.5)
segments(10, ibc_m$m_10_lts[6], 15, ibc_m$m_15_lts[6], col=brew[6], lwd=1.5)
segments(15, ibc_m$m_15_lts[6], 20, ibc_m$m_20_lts[6], col=brew[6], lwd=1.5)
segments(20, ibc_m$m_20_lts[6], 25, ibc_m$m_25_lts[6], col=brew[6], lwd=1.5)
segments(25, ibc_m$m_25_lts[6], 30, ibc_m$m_30_lts[6], col=brew[6], lwd=1.5)
segments(30, ibc_m$m_30_lts[6], 35, ibc_m$m_35_lts[6], col=brew[6], lwd=1.5)
segments(35, ibc_m$m_35_lts[6], 40, ibc_m$m_40_lts[6], col=brew[6], lwd=1.5)
segments(40, ibc_m$m_40_lts[6], 45, ibc_m$m_45_lts[6], col=brew[6], lwd=1.5)
segments(45, ibc_m$m_45_lts[6], 50, ibc_m$m_50_lts[6], col=brew[6], lwd=1.5)
segments(50, ibc_m$m_50_lts[6], 55, ibc_m$m_55_lts[6], col=brew[6], lwd=1.5)
segments(55, ibc_m$m_55_lts[6], 60, ibc_m$m_60_lts[6], col=brew[6], lwd=1.5)
#SVMPO
lines(5, ibc_m$m_05_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(10, ibc_m$m_10_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(15, ibc_m$m_15_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(20, ibc_m$m_20_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(25, ibc_m$m_25_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(30, ibc_m$m_30_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(35, ibc_m$m_35_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(40, ibc_m$m_40_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(45, ibc_m$m_45_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(50, ibc_m$m_50_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(55, ibc_m$m_55_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(60, ibc_m$m_60_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
segments(0, 0, 5, ibc_m$m_05_lts[7], col=brew[7], lwd=1.5)
segments(5,  ibc_m$m_05_lts[7], 10, ibc_m$m_10_lts[7], col=brew[7], lwd=1.5)
segments(10, ibc_m$m_10_lts[7], 15, ibc_m$m_15_lts[7], col=brew[7], lwd=1.5)
segments(15, ibc_m$m_15_lts[7], 20, ibc_m$m_20_lts[7], col=brew[7], lwd=1.5)
segments(20, ibc_m$m_20_lts[7], 25, ibc_m$m_25_lts[7], col=brew[7], lwd=1.5)
segments(25, ibc_m$m_25_lts[7], 30, ibc_m$m_30_lts[7], col=brew[7], lwd=1.5)
segments(30, ibc_m$m_30_lts[7], 35, ibc_m$m_35_lts[7], col=brew[7], lwd=1.5)
segments(35, ibc_m$m_35_lts[7], 40, ibc_m$m_40_lts[7], col=brew[7], lwd=1.5)
segments(40, ibc_m$m_40_lts[7], 45, ibc_m$m_45_lts[7], col=brew[7], lwd=1.5)
segments(45, ibc_m$m_45_lts[7], 50, ibc_m$m_50_lts[7], col=brew[7], lwd=1.5)
segments(50, ibc_m$m_50_lts[7], 55, ibc_m$m_55_lts[7], col=brew[7], lwd=1.5)
segments(55, ibc_m$m_55_lts[7], 60, ibc_m$m_60_lts[7], col=brew[7], lwd=1.5)
#YMPO
lines(5, ibc_m$m_05_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(10, ibc_m$m_10_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(15, ibc_m$m_15_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(20, ibc_m$m_20_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(25, ibc_m$m_25_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(30, ibc_m$m_30_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(35, ibc_m$m_35_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(40, ibc_m$m_40_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(45, ibc_m$m_45_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(50, ibc_m$m_50_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(55, ibc_m$m_55_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(60, ibc_m$m_60_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
segments(0, 0, 5, ibc_m$m_05_lts[8], col=brew[8], lwd=1.5)
segments(5,  ibc_m$m_05_lts[8], 10, ibc_m$m_10_lts[8], col=brew[8], lwd=1.5)
segments(10, ibc_m$m_10_lts[8], 15, ibc_m$m_15_lts[8], col=brew[8], lwd=1.5)
segments(15, ibc_m$m_15_lts[8], 20, ibc_m$m_20_lts[8], col=brew[8], lwd=1.5)
segments(20, ibc_m$m_20_lts[8], 25, ibc_m$m_25_lts[8], col=brew[8], lwd=1.5)
segments(25, ibc_m$m_25_lts[8], 30, ibc_m$m_30_lts[8], col=brew[8], lwd=1.5)
segments(30, ibc_m$m_30_lts[8], 35, ibc_m$m_35_lts[8], col=brew[8], lwd=1.5)
segments(35, ibc_m$m_35_lts[8], 40, ibc_m$m_40_lts[8], col=brew[8], lwd=1.5)
segments(40, ibc_m$m_40_lts[8], 45, ibc_m$m_45_lts[8], col=brew[8], lwd=1.5)
segments(45, ibc_m$m_45_lts[8], 50, ibc_m$m_50_lts[8], col=brew[8], lwd=1.5)
segments(50, ibc_m$m_50_lts[8], 55, ibc_m$m_55_lts[8], col=brew[8], lwd=1.5)
segments(55, ibc_m$m_55_lts[8], 60, ibc_m$m_60_lts[8], col=brew[8], lwd=1.5)
# Create x-axis
axis(side=1, at=c(0,5,10,15,20,25,30,35,40,45,50,55,60), tick=TRUE, cex=1.2)
# Create legend
text(0, 10, "Legend: IBC Cyclist", adj=0, cex=0.8, font=2, col="black")
points(0, 9.6, col=brew[1], cex=0.8, pch=20)
segments(0, 9.6, 5, 9.6, col=brew[1], lty=1, lwd=1.5)
text(7, 9.6, "CYMPO", adj=0, cex=0.8, col="black")
points(0, 9.4, col=brew[3], cex=0.8, pch=20)
segments(0, 9.4, 5, 9.4, col=brew[3], lty=1, lwd=1.5)
text(7, 9.4, "LHMPO", adj=0, cex=0.8, col="black")
points(0, 9.2, col=brew[4], cex=0.8, pch=20)
segments(0, 9.2, 5, 9.2, col=brew[4], lty=1, lwd=1.5)
text(7, 9.2, "MAG", adj=0, cex=0.8, col="black")
points(0, 9.0, col=brew[2], cex=0.8, pch=20)
segments(0, 9.0, 5, 9.0, col=brew[2], lty=1, lwd=1.5)
text(7, 9.0, "MetroPlan", adj=0, cex=0.8, col="black")
points(0, 8.8, col=brew[5], cex=0.8, pch=20)
segments(0, 8.8, 5, 8.8, col=brew[5], lty=1, lwd=1.5)
text(7, 8.8, "PAG", adj=0, cex=0.8, col="black")
points(0, 8.6, col=brew[6], cex=0.8, pch=20)
segments(0, 8.6, 5, 8.6, col=brew[6], lty=1, lwd=1.5)
text(7, 8.6, "SCMPO", adj=0, cex=0.8, col="black")
points(0, 8.4, col=brew[7], cex=0.8, pch=20)
segments(0, 8.4, 5, 8.4, col=brew[7], lty=1, lwd=1.5)
text(7, 8.4, "SVMPO", adj=0, cex=0.8, col="black")
points(0, 8.2, col=brew[8], cex=0.8, pch=20)
segments(0, 8.2, 5, 8.2, col=brew[8], lty=1, lwd=1.5)
text(7, 8.2, "YMPO", adj=0, cex=0.8, col="black")
rect(-1, 8.0, 14, 10.2, border="black", col=NA)
graphics.off()

### EAC Cyclist Type

# Graphic set-up
png(file=paste(dat_dir2, "plot_eac-markets.png", sep="/_analysis/_plots/"), bg="white", width=9, height=9, units='in', res=750)
plot(c(0,60), c(0,10), type="n", xlab="Travel time (minutes)", ylab="Market accessibility", xaxt="n", cex=1.2)
# Plot points and lines
# CYMPO
lines(5, eac_m$m_05_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(10, eac_m$m_10_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(15, eac_m$m_15_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(20, eac_m$m_20_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(25, eac_m$m_25_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(30, eac_m$m_30_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(35, eac_m$m_35_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(40, eac_m$m_40_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(45, eac_m$m_45_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(50, eac_m$m_50_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(55, eac_m$m_55_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(60, eac_m$m_60_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
segments(0, 0, 5, eac_m$m_05_lts[1], col=brew[1], lwd=1.5)
segments(5,  eac_m$m_05_lts[1], 10, eac_m$m_10_lts[1], col=brew[1], lwd=1.5)
segments(10, eac_m$m_10_lts[1], 15, eac_m$m_15_lts[1], col=brew[1], lwd=1.5)
segments(15, eac_m$m_15_lts[1], 20, eac_m$m_20_lts[1], col=brew[1], lwd=1.5)
segments(20, eac_m$m_20_lts[1], 25, eac_m$m_25_lts[1], col=brew[1], lwd=1.5)
segments(25, eac_m$m_25_lts[1], 30, eac_m$m_30_lts[1], col=brew[1], lwd=1.5)
segments(30, eac_m$m_30_lts[1], 35, eac_m$m_35_lts[1], col=brew[1], lwd=1.5)
segments(35, eac_m$m_35_lts[1], 40, eac_m$m_40_lts[1], col=brew[1], lwd=1.5)
segments(40, eac_m$m_40_lts[1], 45, eac_m$m_45_lts[1], col=brew[1], lwd=1.5)
segments(45, eac_m$m_45_lts[1], 50, eac_m$m_50_lts[1], col=brew[1], lwd=1.5)
segments(50, eac_m$m_50_lts[1], 55, eac_m$m_55_lts[1], col=brew[1], lwd=1.5)
segments(55, eac_m$m_55_lts[1], 60, eac_m$m_60_lts[1], col=brew[1], lwd=1.5)
#FMPO
lines(5, eac_m$m_05_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(10, eac_m$m_10_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(15, eac_m$m_15_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(20, eac_m$m_20_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(25, eac_m$m_25_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(30, eac_m$m_30_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(35, eac_m$m_35_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(40, eac_m$m_40_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(45, eac_m$m_45_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(50, eac_m$m_50_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(55, eac_m$m_55_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(60, eac_m$m_60_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
segments(0, 0, 5, eac_m$m_05_lts[2], col=brew[2], lwd=1.5)
segments(5,  eac_m$m_05_lts[2], 10, eac_m$m_10_lts[2], col=brew[2], lwd=1.5)
segments(10, eac_m$m_10_lts[2], 15, eac_m$m_15_lts[2], col=brew[2], lwd=1.5)
segments(15, eac_m$m_15_lts[2], 20, eac_m$m_20_lts[2], col=brew[2], lwd=1.5)
segments(20, eac_m$m_20_lts[2], 25, eac_m$m_25_lts[2], col=brew[2], lwd=1.5)
segments(25, eac_m$m_25_lts[2], 30, eac_m$m_30_lts[2], col=brew[2], lwd=1.5)
segments(30, eac_m$m_30_lts[2], 35, eac_m$m_35_lts[2], col=brew[2], lwd=1.5)
segments(35, eac_m$m_35_lts[2], 40, eac_m$m_40_lts[2], col=brew[2], lwd=1.5)
segments(40, eac_m$m_40_lts[2], 45, eac_m$m_45_lts[2], col=brew[2], lwd=1.5)
segments(45, eac_m$m_45_lts[2], 50, eac_m$m_50_lts[2], col=brew[2], lwd=1.5)
segments(50, eac_m$m_50_lts[2], 55, eac_m$m_55_lts[2], col=brew[2], lwd=1.5)
segments(55, eac_m$m_55_lts[2], 60, eac_m$m_60_lts[2], col=brew[2], lwd=1.5)
#LHMPO
lines(5, eac_m$m_05_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(10, eac_m$m_10_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(15, eac_m$m_15_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(20, eac_m$m_20_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(25, eac_m$m_25_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(30, eac_m$m_30_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(35, eac_m$m_35_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(40, eac_m$m_40_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(45, eac_m$m_45_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(50, eac_m$m_50_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(55, eac_m$m_55_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(60, eac_m$m_60_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
segments(0, 0, 5, eac_m$m_05_lts[3], col=brew[3], lwd=1.5)
segments(5,  eac_m$m_05_lts[3], 10, eac_m$m_10_lts[3], col=brew[3], lwd=1.5)
segments(10, eac_m$m_10_lts[3], 15, eac_m$m_15_lts[3], col=brew[3], lwd=1.5)
segments(15, eac_m$m_15_lts[3], 20, eac_m$m_20_lts[3], col=brew[3], lwd=1.5)
segments(20, eac_m$m_20_lts[3], 25, eac_m$m_25_lts[3], col=brew[3], lwd=1.5)
segments(25, eac_m$m_25_lts[3], 30, eac_m$m_30_lts[3], col=brew[3], lwd=1.5)
segments(30, eac_m$m_30_lts[3], 35, eac_m$m_35_lts[3], col=brew[3], lwd=1.5)
segments(35, eac_m$m_35_lts[3], 40, eac_m$m_40_lts[3], col=brew[3], lwd=1.5)
segments(40, eac_m$m_40_lts[3], 45, eac_m$m_45_lts[3], col=brew[3], lwd=1.5)
segments(45, eac_m$m_45_lts[3], 50, eac_m$m_50_lts[3], col=brew[3], lwd=1.5)
segments(50, eac_m$m_50_lts[3], 55, eac_m$m_55_lts[3], col=brew[3], lwd=1.5)
segments(55, eac_m$m_55_lts[3], 60, eac_m$m_60_lts[3], col=brew[3], lwd=1.5)
#MAG
lines(5, eac_m$m_05_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(10, eac_m$m_10_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(15, eac_m$m_15_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(20, eac_m$m_20_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(25, eac_m$m_25_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(30, eac_m$m_30_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(35, eac_m$m_35_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(40, eac_m$m_40_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(45, eac_m$m_45_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(50, eac_m$m_50_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(55, eac_m$m_55_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(60, eac_m$m_60_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
segments(0, 0, 5, eac_m$m_05_lts[4], col=brew[4], lwd=1.5)
segments(5,  eac_m$m_05_lts[4], 10, eac_m$m_10_lts[4], col=brew[4], lwd=1.5)
segments(10, eac_m$m_10_lts[4], 15, eac_m$m_15_lts[4], col=brew[4], lwd=1.5)
segments(15, eac_m$m_15_lts[4], 20, eac_m$m_20_lts[4], col=brew[4], lwd=1.5)
segments(20, eac_m$m_20_lts[4], 25, eac_m$m_25_lts[4], col=brew[4], lwd=1.5)
segments(25, eac_m$m_25_lts[4], 30, eac_m$m_30_lts[4], col=brew[4], lwd=1.5)
segments(30, eac_m$m_30_lts[4], 35, eac_m$m_35_lts[4], col=brew[4], lwd=1.5)
segments(35, eac_m$m_35_lts[4], 40, eac_m$m_40_lts[4], col=brew[4], lwd=1.5)
segments(40, eac_m$m_40_lts[4], 45, eac_m$m_45_lts[4], col=brew[4], lwd=1.5)
segments(45, eac_m$m_45_lts[4], 50, eac_m$m_50_lts[4], col=brew[4], lwd=1.5)
segments(50, eac_m$m_50_lts[4], 55, eac_m$m_55_lts[4], col=brew[4], lwd=1.5)
segments(55, eac_m$m_55_lts[4], 60, eac_m$m_60_lts[4], col=brew[4], lwd=1.5)
#PAG
lines(5, eac_m$m_05_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(10, eac_m$m_10_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(15, eac_m$m_15_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(20, eac_m$m_20_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(25, eac_m$m_25_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(30, eac_m$m_30_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(35, eac_m$m_35_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(40, eac_m$m_40_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(45, eac_m$m_45_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(50, eac_m$m_50_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(55, eac_m$m_55_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(60, eac_m$m_60_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
segments(0, 0, 5, eac_m$m_05_lts[5], col=brew[5], lwd=1.5)
segments(5,  eac_m$m_05_lts[5], 10, eac_m$m_10_lts[5], col=brew[5], lwd=1.5)
segments(10, eac_m$m_10_lts[5], 15, eac_m$m_15_lts[5], col=brew[5], lwd=1.5)
segments(15, eac_m$m_15_lts[5], 20, eac_m$m_20_lts[5], col=brew[5], lwd=1.5)
segments(20, eac_m$m_20_lts[5], 25, eac_m$m_25_lts[5], col=brew[5], lwd=1.5)
segments(25, eac_m$m_25_lts[5], 30, eac_m$m_30_lts[5], col=brew[5], lwd=1.5)
segments(30, eac_m$m_30_lts[5], 35, eac_m$m_35_lts[5], col=brew[5], lwd=1.5)
segments(35, eac_m$m_35_lts[5], 40, eac_m$m_40_lts[5], col=brew[5], lwd=1.5)
segments(40, eac_m$m_40_lts[5], 45, eac_m$m_45_lts[5], col=brew[5], lwd=1.5)
segments(45, eac_m$m_45_lts[5], 50, eac_m$m_50_lts[5], col=brew[5], lwd=1.5)
segments(50, eac_m$m_50_lts[5], 55, eac_m$m_55_lts[5], col=brew[5], lwd=1.5)
segments(55, eac_m$m_55_lts[5], 60, eac_m$m_60_lts[5], col=brew[5], lwd=1.5)
#SCMPO
lines(5, eac_m$m_05_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(10, eac_m$m_10_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(15, eac_m$m_15_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(20, eac_m$m_20_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(25, eac_m$m_25_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(30, eac_m$m_30_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(35, eac_m$m_35_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(40, eac_m$m_40_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(45, eac_m$m_45_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(50, eac_m$m_50_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(55, eac_m$m_55_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(60, eac_m$m_60_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
segments(0, 0, 5, eac_m$m_05_lts[6], col=brew[6], lwd=1.5)
segments(5,  eac_m$m_05_lts[6], 10, eac_m$m_10_lts[6], col=brew[6], lwd=1.5)
segments(10, eac_m$m_10_lts[6], 15, eac_m$m_15_lts[6], col=brew[6], lwd=1.5)
segments(15, eac_m$m_15_lts[6], 20, eac_m$m_20_lts[6], col=brew[6], lwd=1.5)
segments(20, eac_m$m_20_lts[6], 25, eac_m$m_25_lts[6], col=brew[6], lwd=1.5)
segments(25, eac_m$m_25_lts[6], 30, eac_m$m_30_lts[6], col=brew[6], lwd=1.5)
segments(30, eac_m$m_30_lts[6], 35, eac_m$m_35_lts[6], col=brew[6], lwd=1.5)
segments(35, eac_m$m_35_lts[6], 40, eac_m$m_40_lts[6], col=brew[6], lwd=1.5)
segments(40, eac_m$m_40_lts[6], 45, eac_m$m_45_lts[6], col=brew[6], lwd=1.5)
segments(45, eac_m$m_45_lts[6], 50, eac_m$m_50_lts[6], col=brew[6], lwd=1.5)
segments(50, eac_m$m_50_lts[6], 55, eac_m$m_55_lts[6], col=brew[6], lwd=1.5)
segments(55, eac_m$m_55_lts[6], 60, eac_m$m_60_lts[6], col=brew[6], lwd=1.5)
#SVMPO
lines(5, eac_m$m_05_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(10, eac_m$m_10_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(15, eac_m$m_15_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(20, eac_m$m_20_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(25, eac_m$m_25_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(30, eac_m$m_30_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(35, eac_m$m_35_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(40, eac_m$m_40_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(45, eac_m$m_45_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(50, eac_m$m_50_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(55, eac_m$m_55_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(60, eac_m$m_60_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
segments(0, 0, 5, eac_m$m_05_lts[7], col=brew[7], lwd=1.5)
segments(5,  eac_m$m_05_lts[7], 10, eac_m$m_10_lts[7], col=brew[7], lwd=1.5)
segments(10, eac_m$m_10_lts[7], 15, eac_m$m_15_lts[7], col=brew[7], lwd=1.5)
segments(15, eac_m$m_15_lts[7], 20, eac_m$m_20_lts[7], col=brew[7], lwd=1.5)
segments(20, eac_m$m_20_lts[7], 25, eac_m$m_25_lts[7], col=brew[7], lwd=1.5)
segments(25, eac_m$m_25_lts[7], 30, eac_m$m_30_lts[7], col=brew[7], lwd=1.5)
segments(30, eac_m$m_30_lts[7], 35, eac_m$m_35_lts[7], col=brew[7], lwd=1.5)
segments(35, eac_m$m_35_lts[7], 40, eac_m$m_40_lts[7], col=brew[7], lwd=1.5)
segments(40, eac_m$m_40_lts[7], 45, eac_m$m_45_lts[7], col=brew[7], lwd=1.5)
segments(45, eac_m$m_45_lts[7], 50, eac_m$m_50_lts[7], col=brew[7], lwd=1.5)
segments(50, eac_m$m_50_lts[7], 55, eac_m$m_55_lts[7], col=brew[7], lwd=1.5)
segments(55, eac_m$m_55_lts[7], 60, eac_m$m_60_lts[7], col=brew[7], lwd=1.5)
#YMPO
lines(5, eac_m$m_05_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(10, eac_m$m_10_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(15, eac_m$m_15_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(20, eac_m$m_20_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(25, eac_m$m_25_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(30, eac_m$m_30_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(35, eac_m$m_35_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(40, eac_m$m_40_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(45, eac_m$m_45_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(50, eac_m$m_50_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(55, eac_m$m_55_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(60, eac_m$m_60_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
segments(0, 0, 5, eac_m$m_05_lts[8], col=brew[8], lwd=1.5)
segments(5,  eac_m$m_05_lts[8], 10, eac_m$m_10_lts[8], col=brew[8], lwd=1.5)
segments(10, eac_m$m_10_lts[8], 15, eac_m$m_15_lts[8], col=brew[8], lwd=1.5)
segments(15, eac_m$m_15_lts[8], 20, eac_m$m_20_lts[8], col=brew[8], lwd=1.5)
segments(20, eac_m$m_20_lts[8], 25, eac_m$m_25_lts[8], col=brew[8], lwd=1.5)
segments(25, eac_m$m_25_lts[8], 30, eac_m$m_30_lts[8], col=brew[8], lwd=1.5)
segments(30, eac_m$m_30_lts[8], 35, eac_m$m_35_lts[8], col=brew[8], lwd=1.5)
segments(35, eac_m$m_35_lts[8], 40, eac_m$m_40_lts[8], col=brew[8], lwd=1.5)
segments(40, eac_m$m_40_lts[8], 45, eac_m$m_45_lts[8], col=brew[8], lwd=1.5)
segments(45, eac_m$m_45_lts[8], 50, eac_m$m_50_lts[8], col=brew[8], lwd=1.5)
segments(50, eac_m$m_50_lts[8], 55, eac_m$m_55_lts[8], col=brew[8], lwd=1.5)
segments(55, eac_m$m_55_lts[8], 60, eac_m$m_60_lts[8], col=brew[8], lwd=1.5)
# Create x-axis
axis(side=1, at=c(0,5,10,15,20,25,30,35,40,45,50,55,60), tick=TRUE, cex=1.2)
# Create legend
text(0, 10, "Legend: EAC Cyclist", adj=0, cex=0.8, font=2, col="black")
points(0, 9.6, col=brew[1], cex=0.8, pch=20)
segments(0, 9.6, 5, 9.6, col=brew[1], lty=1, lwd=1.5)
text(7, 9.6, "CYMPO", adj=0, cex=0.8, col="black")
points(0, 9.4, col=brew[3], cex=0.8, pch=20)
segments(0, 9.4, 5, 9.4, col=brew[3], lty=1, lwd=1.5)
text(7, 9.4, "LHMPO", adj=0, cex=0.8, col="black")
points(0, 9.2, col=brew[4], cex=0.8, pch=20)
segments(0, 9.2, 5, 9.2, col=brew[4], lty=1, lwd=1.5)
text(7, 9.2, "MAG", adj=0, cex=0.8, col="black")
points(0, 9.0, col=brew[2], cex=0.8, pch=20)
segments(0, 9.0, 5, 9.0, col=brew[2], lty=1, lwd=1.5)
text(7, 9.0, "MetroPlan", adj=0, cex=0.8, col="black")
points(0, 8.8, col=brew[5], cex=0.8, pch=20)
segments(0, 8.8, 5, 8.8, col=brew[5], lty=1, lwd=1.5)
text(7, 8.8, "PAG", adj=0, cex=0.8, col="black")
points(0, 8.6, col=brew[6], cex=0.8, pch=20)
segments(0, 8.6, 5, 8.6, col=brew[6], lty=1, lwd=1.5)
text(7, 8.6, "SCMPO", adj=0, cex=0.8, col="black")
points(0, 8.4, col=brew[7], cex=0.8, pch=20)
segments(0, 8.4, 5, 8.4, col=brew[7], lty=1, lwd=1.5)
text(7, 8.4, "SVMPO", adj=0, cex=0.8, col="black")
points(0, 8.2, col=brew[8], cex=0.8, pch=20)
segments(0, 8.2, 5, 8.2, col=brew[8], lty=1, lwd=1.5)
text(7, 8.2, "YMPO", adj=0, cex=0.8, col="black")
rect(-1, 8.0, 14, 10.2, border="black", col=NA)
graphics.off()

### SAF Cyclist Type

# Graphic set-up
png(file=paste(dat_dir2, "plot_saf-markets.png", sep="/_analysis/_plots/"), bg="white", width=9, height=9, units='in', res=750)
plot(c(0,60), c(0,10), type="n", xlab="Travel time (minutes)", ylab="Market accessibility", xaxt="n", cex=1.2)
# Plot points and lines
# CYMPO
lines(5, saf_m$m_05_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(10, saf_m$m_10_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(15, saf_m$m_15_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(20, saf_m$m_20_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(25, saf_m$m_25_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(30, saf_m$m_30_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(35, saf_m$m_35_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(40, saf_m$m_40_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(45, saf_m$m_45_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(50, saf_m$m_50_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(55, saf_m$m_55_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(60, saf_m$m_60_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
segments(0, 0, 5, saf_m$m_05_lts[1], col=brew[1], lwd=1.5)
segments(5,  saf_m$m_05_lts[1], 10, saf_m$m_10_lts[1], col=brew[1], lwd=1.5)
segments(10, saf_m$m_10_lts[1], 15, saf_m$m_15_lts[1], col=brew[1], lwd=1.5)
segments(15, saf_m$m_15_lts[1], 20, saf_m$m_20_lts[1], col=brew[1], lwd=1.5)
segments(20, saf_m$m_20_lts[1], 25, saf_m$m_25_lts[1], col=brew[1], lwd=1.5)
segments(25, saf_m$m_25_lts[1], 30, saf_m$m_30_lts[1], col=brew[1], lwd=1.5)
segments(30, saf_m$m_30_lts[1], 35, saf_m$m_35_lts[1], col=brew[1], lwd=1.5)
segments(35, saf_m$m_35_lts[1], 40, saf_m$m_40_lts[1], col=brew[1], lwd=1.5)
segments(40, saf_m$m_40_lts[1], 45, saf_m$m_45_lts[1], col=brew[1], lwd=1.5)
segments(45, saf_m$m_45_lts[1], 50, saf_m$m_50_lts[1], col=brew[1], lwd=1.5)
segments(50, saf_m$m_50_lts[1], 55, saf_m$m_55_lts[1], col=brew[1], lwd=1.5)
segments(55, saf_m$m_55_lts[1], 60, saf_m$m_60_lts[1], col=brew[1], lwd=1.5)
#FMPO
lines(5, saf_m$m_05_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(10, saf_m$m_10_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(15, saf_m$m_15_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(20, saf_m$m_20_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(25, saf_m$m_25_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(30, saf_m$m_30_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(35, saf_m$m_35_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(40, saf_m$m_40_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(45, saf_m$m_45_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(50, saf_m$m_50_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(55, saf_m$m_55_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(60, saf_m$m_60_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
segments(0, 0, 5, saf_m$m_05_lts[2], col=brew[2], lwd=1.5)
segments(5,  saf_m$m_05_lts[2], 10, saf_m$m_10_lts[2], col=brew[2], lwd=1.5)
segments(10, saf_m$m_10_lts[2], 15, saf_m$m_15_lts[2], col=brew[2], lwd=1.5)
segments(15, saf_m$m_15_lts[2], 20, saf_m$m_20_lts[2], col=brew[2], lwd=1.5)
segments(20, saf_m$m_20_lts[2], 25, saf_m$m_25_lts[2], col=brew[2], lwd=1.5)
segments(25, saf_m$m_25_lts[2], 30, saf_m$m_30_lts[2], col=brew[2], lwd=1.5)
segments(30, saf_m$m_30_lts[2], 35, saf_m$m_35_lts[2], col=brew[2], lwd=1.5)
segments(35, saf_m$m_35_lts[2], 40, saf_m$m_40_lts[2], col=brew[2], lwd=1.5)
segments(40, saf_m$m_40_lts[2], 45, saf_m$m_45_lts[2], col=brew[2], lwd=1.5)
segments(45, saf_m$m_45_lts[2], 50, saf_m$m_50_lts[2], col=brew[2], lwd=1.5)
segments(50, saf_m$m_50_lts[2], 55, saf_m$m_55_lts[2], col=brew[2], lwd=1.5)
segments(55, saf_m$m_55_lts[2], 60, saf_m$m_60_lts[2], col=brew[2], lwd=1.5)
#LHMPO
lines(5, saf_m$m_05_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(10, saf_m$m_10_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(15, saf_m$m_15_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(20, saf_m$m_20_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(25, saf_m$m_25_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(30, saf_m$m_30_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(35, saf_m$m_35_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(40, saf_m$m_40_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(45, saf_m$m_45_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(50, saf_m$m_50_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(55, saf_m$m_55_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(60, saf_m$m_60_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
segments(0, 0, 5, saf_m$m_05_lts[3], col=brew[3], lwd=1.5)
segments(5,  saf_m$m_05_lts[3], 10, saf_m$m_10_lts[3], col=brew[3], lwd=1.5)
segments(10, saf_m$m_10_lts[3], 15, saf_m$m_15_lts[3], col=brew[3], lwd=1.5)
segments(15, saf_m$m_15_lts[3], 20, saf_m$m_20_lts[3], col=brew[3], lwd=1.5)
segments(20, saf_m$m_20_lts[3], 25, saf_m$m_25_lts[3], col=brew[3], lwd=1.5)
segments(25, saf_m$m_25_lts[3], 30, saf_m$m_30_lts[3], col=brew[3], lwd=1.5)
segments(30, saf_m$m_30_lts[3], 35, saf_m$m_35_lts[3], col=brew[3], lwd=1.5)
segments(35, saf_m$m_35_lts[3], 40, saf_m$m_40_lts[3], col=brew[3], lwd=1.5)
segments(40, saf_m$m_40_lts[3], 45, saf_m$m_45_lts[3], col=brew[3], lwd=1.5)
segments(45, saf_m$m_45_lts[3], 50, saf_m$m_50_lts[3], col=brew[3], lwd=1.5)
segments(50, saf_m$m_50_lts[3], 55, saf_m$m_55_lts[3], col=brew[3], lwd=1.5)
segments(55, saf_m$m_55_lts[3], 60, saf_m$m_60_lts[3], col=brew[3], lwd=1.5)
#MAG
lines(5, saf_m$m_05_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(10, saf_m$m_10_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(15, saf_m$m_15_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(20, saf_m$m_20_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(25, saf_m$m_25_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(30, saf_m$m_30_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(35, saf_m$m_35_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(40, saf_m$m_40_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(45, saf_m$m_45_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(50, saf_m$m_50_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(55, saf_m$m_55_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(60, saf_m$m_60_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
segments(0, 0, 5, saf_m$m_05_lts[4], col=brew[4], lwd=1.5)
segments(5,  saf_m$m_05_lts[4], 10, saf_m$m_10_lts[4], col=brew[4], lwd=1.5)
segments(10, saf_m$m_10_lts[4], 15, saf_m$m_15_lts[4], col=brew[4], lwd=1.5)
segments(15, saf_m$m_15_lts[4], 20, saf_m$m_20_lts[4], col=brew[4], lwd=1.5)
segments(20, saf_m$m_20_lts[4], 25, saf_m$m_25_lts[4], col=brew[4], lwd=1.5)
segments(25, saf_m$m_25_lts[4], 30, saf_m$m_30_lts[4], col=brew[4], lwd=1.5)
segments(30, saf_m$m_30_lts[4], 35, saf_m$m_35_lts[4], col=brew[4], lwd=1.5)
segments(35, saf_m$m_35_lts[4], 40, saf_m$m_40_lts[4], col=brew[4], lwd=1.5)
segments(40, saf_m$m_40_lts[4], 45, saf_m$m_45_lts[4], col=brew[4], lwd=1.5)
segments(45, saf_m$m_45_lts[4], 50, saf_m$m_50_lts[4], col=brew[4], lwd=1.5)
segments(50, saf_m$m_50_lts[4], 55, saf_m$m_55_lts[4], col=brew[4], lwd=1.5)
segments(55, saf_m$m_55_lts[4], 60, saf_m$m_60_lts[4], col=brew[4], lwd=1.5)
#PAG
lines(5, saf_m$m_05_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(10, saf_m$m_10_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(15, saf_m$m_15_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(20, saf_m$m_20_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(25, saf_m$m_25_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(30, saf_m$m_30_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(35, saf_m$m_35_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(40, saf_m$m_40_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(45, saf_m$m_45_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(50, saf_m$m_50_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(55, saf_m$m_55_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(60, saf_m$m_60_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
segments(0, 0, 5, saf_m$m_05_lts[5], col=brew[5], lwd=1.5)
segments(5,  saf_m$m_05_lts[5], 10, saf_m$m_10_lts[5], col=brew[5], lwd=1.5)
segments(10, saf_m$m_10_lts[5], 15, saf_m$m_15_lts[5], col=brew[5], lwd=1.5)
segments(15, saf_m$m_15_lts[5], 20, saf_m$m_20_lts[5], col=brew[5], lwd=1.5)
segments(20, saf_m$m_20_lts[5], 25, saf_m$m_25_lts[5], col=brew[5], lwd=1.5)
segments(25, saf_m$m_25_lts[5], 30, saf_m$m_30_lts[5], col=brew[5], lwd=1.5)
segments(30, saf_m$m_30_lts[5], 35, saf_m$m_35_lts[5], col=brew[5], lwd=1.5)
segments(35, saf_m$m_35_lts[5], 40, saf_m$m_40_lts[5], col=brew[5], lwd=1.5)
segments(40, saf_m$m_40_lts[5], 45, saf_m$m_45_lts[5], col=brew[5], lwd=1.5)
segments(45, saf_m$m_45_lts[5], 50, saf_m$m_50_lts[5], col=brew[5], lwd=1.5)
segments(50, saf_m$m_50_lts[5], 55, saf_m$m_55_lts[5], col=brew[5], lwd=1.5)
segments(55, saf_m$m_55_lts[5], 60, saf_m$m_60_lts[5], col=brew[5], lwd=1.5)
#SCMPO
lines(5, saf_m$m_05_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(10, saf_m$m_10_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(15, saf_m$m_15_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(20, saf_m$m_20_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(25, saf_m$m_25_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(30, saf_m$m_30_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(35, saf_m$m_35_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(40, saf_m$m_40_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(45, saf_m$m_45_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(50, saf_m$m_50_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(55, saf_m$m_55_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(60, saf_m$m_60_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
segments(0, 0, 5, saf_m$m_05_lts[6], col=brew[6], lwd=1.5)
segments(5,  saf_m$m_05_lts[6], 10, saf_m$m_10_lts[6], col=brew[6], lwd=1.5)
segments(10, saf_m$m_10_lts[6], 15, saf_m$m_15_lts[6], col=brew[6], lwd=1.5)
segments(15, saf_m$m_15_lts[6], 20, saf_m$m_20_lts[6], col=brew[6], lwd=1.5)
segments(20, saf_m$m_20_lts[6], 25, saf_m$m_25_lts[6], col=brew[6], lwd=1.5)
segments(25, saf_m$m_25_lts[6], 30, saf_m$m_30_lts[6], col=brew[6], lwd=1.5)
segments(30, saf_m$m_30_lts[6], 35, saf_m$m_35_lts[6], col=brew[6], lwd=1.5)
segments(35, saf_m$m_35_lts[6], 40, saf_m$m_40_lts[6], col=brew[6], lwd=1.5)
segments(40, saf_m$m_40_lts[6], 45, saf_m$m_45_lts[6], col=brew[6], lwd=1.5)
segments(45, saf_m$m_45_lts[6], 50, saf_m$m_50_lts[6], col=brew[6], lwd=1.5)
segments(50, saf_m$m_50_lts[6], 55, saf_m$m_55_lts[6], col=brew[6], lwd=1.5)
segments(55, saf_m$m_55_lts[6], 60, saf_m$m_60_lts[6], col=brew[6], lwd=1.5)
#SVMPO
lines(5, saf_m$m_05_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(10, saf_m$m_10_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(15, saf_m$m_15_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(20, saf_m$m_20_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(25, saf_m$m_25_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(30, saf_m$m_30_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(35, saf_m$m_35_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(40, saf_m$m_40_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(45, saf_m$m_45_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(50, saf_m$m_50_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(55, saf_m$m_55_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(60, saf_m$m_60_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
segments(0, 0, 5, saf_m$m_05_lts[7], col=brew[7], lwd=1.5)
segments(5,  saf_m$m_05_lts[7], 10, saf_m$m_10_lts[7], col=brew[7], lwd=1.5)
segments(10, saf_m$m_10_lts[7], 15, saf_m$m_15_lts[7], col=brew[7], lwd=1.5)
segments(15, saf_m$m_15_lts[7], 20, saf_m$m_20_lts[7], col=brew[7], lwd=1.5)
segments(20, saf_m$m_20_lts[7], 25, saf_m$m_25_lts[7], col=brew[7], lwd=1.5)
segments(25, saf_m$m_25_lts[7], 30, saf_m$m_30_lts[7], col=brew[7], lwd=1.5)
segments(30, saf_m$m_30_lts[7], 35, saf_m$m_35_lts[7], col=brew[7], lwd=1.5)
segments(35, saf_m$m_35_lts[7], 40, saf_m$m_40_lts[7], col=brew[7], lwd=1.5)
segments(40, saf_m$m_40_lts[7], 45, saf_m$m_45_lts[7], col=brew[7], lwd=1.5)
segments(45, saf_m$m_45_lts[7], 50, saf_m$m_50_lts[7], col=brew[7], lwd=1.5)
segments(50, saf_m$m_50_lts[7], 55, saf_m$m_55_lts[7], col=brew[7], lwd=1.5)
segments(55, saf_m$m_55_lts[7], 60, saf_m$m_60_lts[7], col=brew[7], lwd=1.5)
#YMPO
lines(5, saf_m$m_05_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(10, saf_m$m_10_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(15, saf_m$m_15_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(20, saf_m$m_20_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(25, saf_m$m_25_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(30, saf_m$m_30_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(35, saf_m$m_35_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(40, saf_m$m_40_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(45, saf_m$m_45_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(50, saf_m$m_50_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(55, saf_m$m_55_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(60, saf_m$m_60_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
segments(0, 0, 5, saf_m$m_05_lts[8], col=brew[8], lwd=1.5)
segments(5,  saf_m$m_05_lts[8], 10, saf_m$m_10_lts[8], col=brew[8], lwd=1.5)
segments(10, saf_m$m_10_lts[8], 15, saf_m$m_15_lts[8], col=brew[8], lwd=1.5)
segments(15, saf_m$m_15_lts[8], 20, saf_m$m_20_lts[8], col=brew[8], lwd=1.5)
segments(20, saf_m$m_20_lts[8], 25, saf_m$m_25_lts[8], col=brew[8], lwd=1.5)
segments(25, saf_m$m_25_lts[8], 30, saf_m$m_30_lts[8], col=brew[8], lwd=1.5)
segments(30, saf_m$m_30_lts[8], 35, saf_m$m_35_lts[8], col=brew[8], lwd=1.5)
segments(35, saf_m$m_35_lts[8], 40, saf_m$m_40_lts[8], col=brew[8], lwd=1.5)
segments(40, saf_m$m_40_lts[8], 45, saf_m$m_45_lts[8], col=brew[8], lwd=1.5)
segments(45, saf_m$m_45_lts[8], 50, saf_m$m_50_lts[8], col=brew[8], lwd=1.5)
segments(50, saf_m$m_50_lts[8], 55, saf_m$m_55_lts[8], col=brew[8], lwd=1.5)
segments(55, saf_m$m_55_lts[8], 60, saf_m$m_60_lts[8], col=brew[8], lwd=1.5)
# Create x-axis
axis(side=1, at=c(0,5,10,15,20,25,30,35,40,45,50,55,60), tick=TRUE, cex=1.2)
# Create legend
text(0, 10, "Legend: SAF Cyclist", adj=0, cex=0.8, font=2, col="black")
points(0, 9.6, col=brew[1], cex=0.8, pch=20)
segments(0, 9.6, 5, 9.6, col=brew[1], lty=1, lwd=1.5)
text(7, 9.6, "CYMPO", adj=0, cex=0.8, col="black")
points(0, 9.4, col=brew[3], cex=0.8, pch=20)
segments(0, 9.4, 5, 9.4, col=brew[3], lty=1, lwd=1.5)
text(7, 9.4, "LHMPO", adj=0, cex=0.8, col="black")
points(0, 9.2, col=brew[4], cex=0.8, pch=20)
segments(0, 9.2, 5, 9.2, col=brew[4], lty=1, lwd=1.5)
text(7, 9.2, "MAG", adj=0, cex=0.8, col="black")
points(0, 9.0, col=brew[2], cex=0.8, pch=20)
segments(0, 9.0, 5, 9.0, col=brew[2], lty=1, lwd=1.5)
text(7, 9.0, "MetroPlan", adj=0, cex=0.8, col="black")
points(0, 8.8, col=brew[5], cex=0.8, pch=20)
segments(0, 8.8, 5, 8.8, col=brew[5], lty=1, lwd=1.5)
text(7, 8.8, "PAG", adj=0, cex=0.8, col="black")
points(0, 8.6, col=brew[6], cex=0.8, pch=20)
segments(0, 8.6, 5, 8.6, col=brew[6], lty=1, lwd=1.5)
text(7, 8.6, "SCMPO", adj=0, cex=0.8, col="black")
points(0, 8.4, col=brew[7], cex=0.8, pch=20)
segments(0, 8.4, 5, 8.4, col=brew[7], lty=1, lwd=1.5)
text(7, 8.4, "SVMPO", adj=0, cex=0.8, col="black")
points(0, 8.2, col=brew[8], cex=0.8, pch=20)
segments(0, 8.2, 5, 8.2, col=brew[8], lty=1, lwd=1.5)
text(7, 8.2, "YMPO", adj=0, cex=0.8, col="black")
rect(-1, 8.0, 14, 10.2, border="black", col=NA)
graphics.off()

####################
### SCHOOLS
####################

# Summarize school access by cyclist type
ibc_s <- dat_ibc_s %>% group_by(mpo) %>%
  summarise(s_05=mean(s_00_05,na.rm=T),s_10=mean(s_05_10,na.rm=T),s_15=mean(s_10_15,na.rm=T),s_20=mean(s_15_20,na.rm=T),
            s_25=mean(s_20_25,na.rm=T),s_30=mean(s_25_30,na.rm=T),s_35=mean(s_30_35,na.rm=T),s_40=mean(s_35_40,na.rm=T),
            s_45=mean(s_40_45,na.rm=T),s_50=mean(s_45_50,na.rm=T),s_55=mean(s_50_55,na.rm=T),s_60=mean(s_55_60,na.rm=T),
            s_05_lts=mean(s_00_05_lt,na.rm=T),s_10_lts=mean(s_05_10_lt,na.rm=T),s_15_lts=mean(s_10_15_lt,na.rm=T),s_20_lts=mean(s_15_20_lt,na.rm=T),
            s_25_lts=mean(s_20_25_lt,na.rm=T),s_30_lts=mean(s_25_30_lt,na.rm=T),s_35_lts=mean(s_30_35_lt,na.rm=T),s_40_lts=mean(s_35_40_lt,na.rm=T),
            s_45_lts=mean(s_40_45_lt,na.rm=T),s_50_lts=mean(s_45_50_lt,na.rm=T),s_55_lts=mean(s_50_55_lt,na.rm=T),s_60_lts=mean(s_55_60_lt,na.rm=T))
eac_s <- dat_eac_s %>% group_by(mpo) %>%
  summarise(s_05=mean(s_00_05,na.rm=T),s_10=mean(s_05_10,na.rm=T),s_15=mean(s_10_15,na.rm=T),s_20=mean(s_15_20,na.rm=T),
            s_25=mean(s_20_25,na.rm=T),s_30=mean(s_25_30,na.rm=T),s_35=mean(s_30_35,na.rm=T),s_40=mean(s_35_40,na.rm=T),
            s_45=mean(s_40_45,na.rm=T),s_50=mean(s_45_50,na.rm=T),s_55=mean(s_50_55,na.rm=T),s_60=mean(s_55_60,na.rm=T),
            s_05_lts=mean(s_00_05_lt,na.rm=T),s_10_lts=mean(s_05_10_lt,na.rm=T),s_15_lts=mean(s_10_15_lt,na.rm=T),s_20_lts=mean(s_15_20_lt,na.rm=T),
            s_25_lts=mean(s_20_25_lt,na.rm=T),s_30_lts=mean(s_25_30_lt,na.rm=T),s_35_lts=mean(s_30_35_lt,na.rm=T),s_40_lts=mean(s_35_40_lt,na.rm=T),
            s_45_lts=mean(s_40_45_lt,na.rm=T),s_50_lts=mean(s_45_50_lt,na.rm=T),s_55_lts=mean(s_50_55_lt,na.rm=T),s_60_lts=mean(s_55_60_lt,na.rm=T))
saf_s <- dat_saf_s %>% group_by(mpo) %>%
  summarise(s_05=mean(s_00_05,na.rm=T),s_10=mean(s_05_10,na.rm=T),s_15=mean(s_10_15,na.rm=T),s_20=mean(s_15_20,na.rm=T),
            s_25=mean(s_20_25,na.rm=T),s_30=mean(s_25_30,na.rm=T),s_35=mean(s_30_35,na.rm=T),s_40=mean(s_35_40,na.rm=T),
            s_45=mean(s_40_45,na.rm=T),s_50=mean(s_45_50,na.rm=T),s_55=mean(s_50_55,na.rm=T),s_60=mean(s_55_60,na.rm=T),
            s_05_lts=mean(s_00_05_lt,na.rm=T),s_10_lts=mean(s_05_10_lt,na.rm=T),s_15_lts=mean(s_10_15_lt,na.rm=T),s_20_lts=mean(s_15_20_lt,na.rm=T),
            s_25_lts=mean(s_20_25_lt,na.rm=T),s_30_lts=mean(s_25_30_lt,na.rm=T),s_35_lts=mean(s_30_35_lt,na.rm=T),s_40_lts=mean(s_35_40_lt,na.rm=T),
            s_45_lts=mean(s_40_45_lt,na.rm=T),s_50_lts=mean(s_45_50_lt,na.rm=T),s_55_lts=mean(s_50_55_lt,na.rm=T),s_60_lts=mean(s_55_60_lt,na.rm=T))

### IBC Cyclist Type

# Graphic set-up
png(file=paste(dat_dir2, "plot_ibc-schools.png", sep="/_analysis/_plots/"), bg="white", width=9, height=9, units='in', res=750)
plot(c(0,60), c(0,25), type="n", xlab="Travel time (minutes)", ylab="School accessibility", xaxt="n", cex=1.2)
# Plot points and lines
# CYMPO
lines(5, ibc_s$s_05_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(10, ibc_s$s_10_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(15, ibc_s$s_15_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(20, ibc_s$s_20_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(25, ibc_s$s_25_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(30, ibc_s$s_30_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(35, ibc_s$s_35_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(40, ibc_s$s_40_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(45, ibc_s$s_45_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(50, ibc_s$s_50_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(55, ibc_s$s_55_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(60, ibc_s$s_60_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
segments(0, 0, 5, ibc_s$s_05_lts[1], col=brew[1], lwd=1.5)
segments(5,  ibc_s$s_05_lts[1], 10, ibc_s$s_10_lts[1], col=brew[1], lwd=1.5)
segments(10, ibc_s$s_10_lts[1], 15, ibc_s$s_15_lts[1], col=brew[1], lwd=1.5)
segments(15, ibc_s$s_15_lts[1], 20, ibc_s$s_20_lts[1], col=brew[1], lwd=1.5)
segments(20, ibc_s$s_20_lts[1], 25, ibc_s$s_25_lts[1], col=brew[1], lwd=1.5)
segments(25, ibc_s$s_25_lts[1], 30, ibc_s$s_30_lts[1], col=brew[1], lwd=1.5)
segments(30, ibc_s$s_30_lts[1], 35, ibc_s$s_35_lts[1], col=brew[1], lwd=1.5)
segments(35, ibc_s$s_35_lts[1], 40, ibc_s$s_40_lts[1], col=brew[1], lwd=1.5)
segments(40, ibc_s$s_40_lts[1], 45, ibc_s$s_45_lts[1], col=brew[1], lwd=1.5)
segments(45, ibc_s$s_45_lts[1], 50, ibc_s$s_50_lts[1], col=brew[1], lwd=1.5)
segments(50, ibc_s$s_50_lts[1], 55, ibc_s$s_55_lts[1], col=brew[1], lwd=1.5)
segments(55, ibc_s$s_55_lts[1], 60, ibc_s$s_60_lts[1], col=brew[1], lwd=1.5)
#FMPO
lines(5, ibc_s$s_05_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(10, ibc_s$s_10_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(15, ibc_s$s_15_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(20, ibc_s$s_20_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(25, ibc_s$s_25_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(30, ibc_s$s_30_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(35, ibc_s$s_35_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(40, ibc_s$s_40_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(45, ibc_s$s_45_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(50, ibc_s$s_50_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(55, ibc_s$s_55_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(60, ibc_s$s_60_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
segments(0, 0, 5, ibc_s$s_05_lts[2], col=brew[2], lwd=1.5)
segments(5,  ibc_s$s_05_lts[2], 10, ibc_s$s_10_lts[2], col=brew[2], lwd=1.5)
segments(10, ibc_s$s_10_lts[2], 15, ibc_s$s_15_lts[2], col=brew[2], lwd=1.5)
segments(15, ibc_s$s_15_lts[2], 20, ibc_s$s_20_lts[2], col=brew[2], lwd=1.5)
segments(20, ibc_s$s_20_lts[2], 25, ibc_s$s_25_lts[2], col=brew[2], lwd=1.5)
segments(25, ibc_s$s_25_lts[2], 30, ibc_s$s_30_lts[2], col=brew[2], lwd=1.5)
segments(30, ibc_s$s_30_lts[2], 35, ibc_s$s_35_lts[2], col=brew[2], lwd=1.5)
segments(35, ibc_s$s_35_lts[2], 40, ibc_s$s_40_lts[2], col=brew[2], lwd=1.5)
segments(40, ibc_s$s_40_lts[2], 45, ibc_s$s_45_lts[2], col=brew[2], lwd=1.5)
segments(45, ibc_s$s_45_lts[2], 50, ibc_s$s_50_lts[2], col=brew[2], lwd=1.5)
segments(50, ibc_s$s_50_lts[2], 55, ibc_s$s_55_lts[2], col=brew[2], lwd=1.5)
segments(55, ibc_s$s_55_lts[2], 60, ibc_s$s_60_lts[2], col=brew[2], lwd=1.5)
#LHMPO
lines(5, ibc_s$s_05_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(10, ibc_s$s_10_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(15, ibc_s$s_15_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(20, ibc_s$s_20_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(25, ibc_s$s_25_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(30, ibc_s$s_30_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(35, ibc_s$s_35_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(40, ibc_s$s_40_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(45, ibc_s$s_45_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(50, ibc_s$s_50_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(55, ibc_s$s_55_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(60, ibc_s$s_60_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
segments(0, 0, 5, ibc_s$s_05_lts[3], col=brew[3], lwd=1.5)
segments(5,  ibc_s$s_05_lts[3], 10, ibc_s$s_10_lts[3], col=brew[3], lwd=1.5)
segments(10, ibc_s$s_10_lts[3], 15, ibc_s$s_15_lts[3], col=brew[3], lwd=1.5)
segments(15, ibc_s$s_15_lts[3], 20, ibc_s$s_20_lts[3], col=brew[3], lwd=1.5)
segments(20, ibc_s$s_20_lts[3], 25, ibc_s$s_25_lts[3], col=brew[3], lwd=1.5)
segments(25, ibc_s$s_25_lts[3], 30, ibc_s$s_30_lts[3], col=brew[3], lwd=1.5)
segments(30, ibc_s$s_30_lts[3], 35, ibc_s$s_35_lts[3], col=brew[3], lwd=1.5)
segments(35, ibc_s$s_35_lts[3], 40, ibc_s$s_40_lts[3], col=brew[3], lwd=1.5)
segments(40, ibc_s$s_40_lts[3], 45, ibc_s$s_45_lts[3], col=brew[3], lwd=1.5)
segments(45, ibc_s$s_45_lts[3], 50, ibc_s$s_50_lts[3], col=brew[3], lwd=1.5)
segments(50, ibc_s$s_50_lts[3], 55, ibc_s$s_55_lts[3], col=brew[3], lwd=1.5)
segments(55, ibc_s$s_55_lts[3], 60, ibc_s$s_60_lts[3], col=brew[3], lwd=1.5)
#MAG
lines(5, ibc_s$s_05_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(10, ibc_s$s_10_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(15, ibc_s$s_15_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(20, ibc_s$s_20_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(25, ibc_s$s_25_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(30, ibc_s$s_30_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(35, ibc_s$s_35_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(40, ibc_s$s_40_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(45, ibc_s$s_45_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(50, ibc_s$s_50_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(55, ibc_s$s_55_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(60, ibc_s$s_60_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
segments(0, 0, 5, ibc_s$s_05_lts[4], col=brew[4], lwd=1.5)
segments(5,  ibc_s$s_05_lts[4], 10, ibc_s$s_10_lts[4], col=brew[4], lwd=1.5)
segments(10, ibc_s$s_10_lts[4], 15, ibc_s$s_15_lts[4], col=brew[4], lwd=1.5)
segments(15, ibc_s$s_15_lts[4], 20, ibc_s$s_20_lts[4], col=brew[4], lwd=1.5)
segments(20, ibc_s$s_20_lts[4], 25, ibc_s$s_25_lts[4], col=brew[4], lwd=1.5)
segments(25, ibc_s$s_25_lts[4], 30, ibc_s$s_30_lts[4], col=brew[4], lwd=1.5)
segments(30, ibc_s$s_30_lts[4], 35, ibc_s$s_35_lts[4], col=brew[4], lwd=1.5)
segments(35, ibc_s$s_35_lts[4], 40, ibc_s$s_40_lts[4], col=brew[4], lwd=1.5)
segments(40, ibc_s$s_40_lts[4], 45, ibc_s$s_45_lts[4], col=brew[4], lwd=1.5)
segments(45, ibc_s$s_45_lts[4], 50, ibc_s$s_50_lts[4], col=brew[4], lwd=1.5)
segments(50, ibc_s$s_50_lts[4], 55, ibc_s$s_55_lts[4], col=brew[4], lwd=1.5)
segments(55, ibc_s$s_55_lts[4], 60, ibc_s$s_60_lts[4], col=brew[4], lwd=1.5)
#PAG
lines(5, ibc_s$s_05_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(10, ibc_s$s_10_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(15, ibc_s$s_15_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(20, ibc_s$s_20_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(25, ibc_s$s_25_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(30, ibc_s$s_30_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(35, ibc_s$s_35_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(40, ibc_s$s_40_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(45, ibc_s$s_45_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(50, ibc_s$s_50_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(55, ibc_s$s_55_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(60, ibc_s$s_60_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
segments(0, 0, 5, ibc_s$s_05_lts[5], col=brew[5], lwd=1.5)
segments(5,  ibc_s$s_05_lts[5], 10, ibc_s$s_10_lts[5], col=brew[5], lwd=1.5)
segments(10, ibc_s$s_10_lts[5], 15, ibc_s$s_15_lts[5], col=brew[5], lwd=1.5)
segments(15, ibc_s$s_15_lts[5], 20, ibc_s$s_20_lts[5], col=brew[5], lwd=1.5)
segments(20, ibc_s$s_20_lts[5], 25, ibc_s$s_25_lts[5], col=brew[5], lwd=1.5)
segments(25, ibc_s$s_25_lts[5], 30, ibc_s$s_30_lts[5], col=brew[5], lwd=1.5)
segments(30, ibc_s$s_30_lts[5], 35, ibc_s$s_35_lts[5], col=brew[5], lwd=1.5)
segments(35, ibc_s$s_35_lts[5], 40, ibc_s$s_40_lts[5], col=brew[5], lwd=1.5)
segments(40, ibc_s$s_40_lts[5], 45, ibc_s$s_45_lts[5], col=brew[5], lwd=1.5)
segments(45, ibc_s$s_45_lts[5], 50, ibc_s$s_50_lts[5], col=brew[5], lwd=1.5)
segments(50, ibc_s$s_50_lts[5], 55, ibc_s$s_55_lts[5], col=brew[5], lwd=1.5)
segments(55, ibc_s$s_55_lts[5], 60, ibc_s$s_60_lts[5], col=brew[5], lwd=1.5)
#SCMPO
lines(5, ibc_s$s_05_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(10, ibc_s$s_10_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(15, ibc_s$s_15_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(20, ibc_s$s_20_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(25, ibc_s$s_25_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(30, ibc_s$s_30_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(35, ibc_s$s_35_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(40, ibc_s$s_40_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(45, ibc_s$s_45_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(50, ibc_s$s_50_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(55, ibc_s$s_55_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(60, ibc_s$s_60_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
segments(0, 0, 5, ibc_s$s_05_lts[6], col=brew[6], lwd=1.5)
segments(5,  ibc_s$s_05_lts[6], 10, ibc_s$s_10_lts[6], col=brew[6], lwd=1.5)
segments(10, ibc_s$s_10_lts[6], 15, ibc_s$s_15_lts[6], col=brew[6], lwd=1.5)
segments(15, ibc_s$s_15_lts[6], 20, ibc_s$s_20_lts[6], col=brew[6], lwd=1.5)
segments(20, ibc_s$s_20_lts[6], 25, ibc_s$s_25_lts[6], col=brew[6], lwd=1.5)
segments(25, ibc_s$s_25_lts[6], 30, ibc_s$s_30_lts[6], col=brew[6], lwd=1.5)
segments(30, ibc_s$s_30_lts[6], 35, ibc_s$s_35_lts[6], col=brew[6], lwd=1.5)
segments(35, ibc_s$s_35_lts[6], 40, ibc_s$s_40_lts[6], col=brew[6], lwd=1.5)
segments(40, ibc_s$s_40_lts[6], 45, ibc_s$s_45_lts[6], col=brew[6], lwd=1.5)
segments(45, ibc_s$s_45_lts[6], 50, ibc_s$s_50_lts[6], col=brew[6], lwd=1.5)
segments(50, ibc_s$s_50_lts[6], 55, ibc_s$s_55_lts[6], col=brew[6], lwd=1.5)
segments(55, ibc_s$s_55_lts[6], 60, ibc_s$s_60_lts[6], col=brew[6], lwd=1.5)
#SVMPO
lines(5, ibc_s$s_05_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(10, ibc_s$s_10_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(15, ibc_s$s_15_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(20, ibc_s$s_20_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(25, ibc_s$s_25_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(30, ibc_s$s_30_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(35, ibc_s$s_35_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(40, ibc_s$s_40_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(45, ibc_s$s_45_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(50, ibc_s$s_50_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(55, ibc_s$s_55_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(60, ibc_s$s_60_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
segments(0, 0, 5, ibc_s$s_05_lts[7], col=brew[7], lwd=1.5)
segments(5,  ibc_s$s_05_lts[7], 10, ibc_s$s_10_lts[7], col=brew[7], lwd=1.5)
segments(10, ibc_s$s_10_lts[7], 15, ibc_s$s_15_lts[7], col=brew[7], lwd=1.5)
segments(15, ibc_s$s_15_lts[7], 20, ibc_s$s_20_lts[7], col=brew[7], lwd=1.5)
segments(20, ibc_s$s_20_lts[7], 25, ibc_s$s_25_lts[7], col=brew[7], lwd=1.5)
segments(25, ibc_s$s_25_lts[7], 30, ibc_s$s_30_lts[7], col=brew[7], lwd=1.5)
segments(30, ibc_s$s_30_lts[7], 35, ibc_s$s_35_lts[7], col=brew[7], lwd=1.5)
segments(35, ibc_s$s_35_lts[7], 40, ibc_s$s_40_lts[7], col=brew[7], lwd=1.5)
segments(40, ibc_s$s_40_lts[7], 45, ibc_s$s_45_lts[7], col=brew[7], lwd=1.5)
segments(45, ibc_s$s_45_lts[7], 50, ibc_s$s_50_lts[7], col=brew[7], lwd=1.5)
segments(50, ibc_s$s_50_lts[7], 55, ibc_s$s_55_lts[7], col=brew[7], lwd=1.5)
segments(55, ibc_s$s_55_lts[7], 60, ibc_s$s_60_lts[7], col=brew[7], lwd=1.5)
#YMPO
lines(5, ibc_s$s_05_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(10, ibc_s$s_10_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(15, ibc_s$s_15_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(20, ibc_s$s_20_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(25, ibc_s$s_25_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(30, ibc_s$s_30_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(35, ibc_s$s_35_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(40, ibc_s$s_40_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(45, ibc_s$s_45_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(50, ibc_s$s_50_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(55, ibc_s$s_55_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(60, ibc_s$s_60_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
segments(0, 0, 5, ibc_s$s_05_lts[8], col=brew[8], lwd=1.5)
segments(5,  ibc_s$s_05_lts[8], 10, ibc_s$s_10_lts[8], col=brew[8], lwd=1.5)
segments(10, ibc_s$s_10_lts[8], 15, ibc_s$s_15_lts[8], col=brew[8], lwd=1.5)
segments(15, ibc_s$s_15_lts[8], 20, ibc_s$s_20_lts[8], col=brew[8], lwd=1.5)
segments(20, ibc_s$s_20_lts[8], 25, ibc_s$s_25_lts[8], col=brew[8], lwd=1.5)
segments(25, ibc_s$s_25_lts[8], 30, ibc_s$s_30_lts[8], col=brew[8], lwd=1.5)
segments(30, ibc_s$s_30_lts[8], 35, ibc_s$s_35_lts[8], col=brew[8], lwd=1.5)
segments(35, ibc_s$s_35_lts[8], 40, ibc_s$s_40_lts[8], col=brew[8], lwd=1.5)
segments(40, ibc_s$s_40_lts[8], 45, ibc_s$s_45_lts[8], col=brew[8], lwd=1.5)
segments(45, ibc_s$s_45_lts[8], 50, ibc_s$s_50_lts[8], col=brew[8], lwd=1.5)
segments(50, ibc_s$s_50_lts[8], 55, ibc_s$s_55_lts[8], col=brew[8], lwd=1.5)
segments(55, ibc_s$s_55_lts[8], 60, ibc_s$s_60_lts[8], col=brew[8], lwd=1.5)
# Create x-axis
axis(side=1, at=c(0,5,10,15,20,25,30,35,40,45,50,55,60), tick=TRUE, cex=1.2)
# Create legend
text(0, 25, "Legend: IBC Cyclist", adj=0, cex=0.8, font=2, col="black")
points(0, 24, col=brew[1], cex=0.8, pch=20)
segments(0, 24, 5, 24, col=brew[1], lty=1, lwd=1.5)
text(7, 24, "CYMPO", adj=0, cex=0.8, col="black")
points(0, 23.5, col=brew[3], cex=0.8, pch=20)
segments(0, 23.5, 5, 23.5, col=brew[3], lty=1, lwd=1.5)
text(7, 23.5, "LHMPO", adj=0, cex=0.8, col="black")
points(0, 23, col=brew[4], cex=0.8, pch=20)
segments(0, 23, 5, 23, col=brew[4], lty=1, lwd=1.5)
text(7, 23, "MAG", adj=0, cex=0.8, col="black")
points(0, 22.5, col=brew[2], cex=0.8, pch=20)
segments(0, 22.5, 5, 22.5, col=brew[2], lty=1, lwd=1.5)
text(7, 22.5, "MetroPlan", adj=0, cex=0.8, col="black")
points(0, 22, col=brew[5], cex=0.8, pch=20)
segments(0, 22, 5, 22, col=brew[5], lty=1, lwd=1.5)
text(7, 22, "PAG", adj=0, cex=0.8, col="black")
points(0, 21.5, col=brew[6], cex=0.8, pch=20)
segments(0, 21.5, 5, 21.5, col=brew[6], lty=1, lwd=1.5)
text(7, 21.5, "SCMPO", adj=0, cex=0.8, col="black")
points(0, 21, col=brew[7], cex=0.8, pch=20)
segments(0, 21, 5, 21, col=brew[7], lty=1, lwd=1.5)
text(7, 21, "SVMPO", adj=0, cex=0.8, col="black")
points(0, 20.5, col=brew[8], cex=0.8, pch=20)
segments(0, 20.5, 5, 20.5, col=brew[8], lty=1, lwd=1.5)
text(7, 20.5, "YMPO", adj=0, cex=0.8, col="black")
rect(-1, 20, 14, 25.5, border="black", col=NA)
graphics.off()

### EAC Cyclist Type

# Graphic set-up
png(file=paste(dat_dir2, "plot_eac-schools.png", sep="/_analysis/_plots/"), bg="white", width=9, height=9, units='in', res=750)
plot(c(0,60), c(0,25), type="n", xlab="Travel time (minutes)", ylab="School accessibility", xaxt="n", cex=1.2)
# Plot points and lines
# CYMPO
lines(5, eac_s$s_05_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(10, eac_s$s_10_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(15, eac_s$s_15_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(20, eac_s$s_20_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(25, eac_s$s_25_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(30, eac_s$s_30_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(35, eac_s$s_35_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(40, eac_s$s_40_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(45, eac_s$s_45_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(50, eac_s$s_50_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(55, eac_s$s_55_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(60, eac_s$s_60_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
segments(0, 0, 5, eac_s$s_05_lts[1], col=brew[1], lwd=1.5)
segments(5,  eac_s$s_05_lts[1], 10, eac_s$s_10_lts[1], col=brew[1], lwd=1.5)
segments(10, eac_s$s_10_lts[1], 15, eac_s$s_15_lts[1], col=brew[1], lwd=1.5)
segments(15, eac_s$s_15_lts[1], 20, eac_s$s_20_lts[1], col=brew[1], lwd=1.5)
segments(20, eac_s$s_20_lts[1], 25, eac_s$s_25_lts[1], col=brew[1], lwd=1.5)
segments(25, eac_s$s_25_lts[1], 30, eac_s$s_30_lts[1], col=brew[1], lwd=1.5)
segments(30, eac_s$s_30_lts[1], 35, eac_s$s_35_lts[1], col=brew[1], lwd=1.5)
segments(35, eac_s$s_35_lts[1], 40, eac_s$s_40_lts[1], col=brew[1], lwd=1.5)
segments(40, eac_s$s_40_lts[1], 45, eac_s$s_45_lts[1], col=brew[1], lwd=1.5)
segments(45, eac_s$s_45_lts[1], 50, eac_s$s_50_lts[1], col=brew[1], lwd=1.5)
segments(50, eac_s$s_50_lts[1], 55, eac_s$s_55_lts[1], col=brew[1], lwd=1.5)
segments(55, eac_s$s_55_lts[1], 60, eac_s$s_60_lts[1], col=brew[1], lwd=1.5)
#FMPO
lines(5, eac_s$s_05_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(10, eac_s$s_10_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(15, eac_s$s_15_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(20, eac_s$s_20_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(25, eac_s$s_25_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(30, eac_s$s_30_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(35, eac_s$s_35_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(40, eac_s$s_40_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(45, eac_s$s_45_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(50, eac_s$s_50_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(55, eac_s$s_55_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(60, eac_s$s_60_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
segments(0, 0, 5, eac_s$s_05_lts[2], col=brew[2], lwd=1.5)
segments(5,  eac_s$s_05_lts[2], 10, eac_s$s_10_lts[2], col=brew[2], lwd=1.5)
segments(10, eac_s$s_10_lts[2], 15, eac_s$s_15_lts[2], col=brew[2], lwd=1.5)
segments(15, eac_s$s_15_lts[2], 20, eac_s$s_20_lts[2], col=brew[2], lwd=1.5)
segments(20, eac_s$s_20_lts[2], 25, eac_s$s_25_lts[2], col=brew[2], lwd=1.5)
segments(25, eac_s$s_25_lts[2], 30, eac_s$s_30_lts[2], col=brew[2], lwd=1.5)
segments(30, eac_s$s_30_lts[2], 35, eac_s$s_35_lts[2], col=brew[2], lwd=1.5)
segments(35, eac_s$s_35_lts[2], 40, eac_s$s_40_lts[2], col=brew[2], lwd=1.5)
segments(40, eac_s$s_40_lts[2], 45, eac_s$s_45_lts[2], col=brew[2], lwd=1.5)
segments(45, eac_s$s_45_lts[2], 50, eac_s$s_50_lts[2], col=brew[2], lwd=1.5)
segments(50, eac_s$s_50_lts[2], 55, eac_s$s_55_lts[2], col=brew[2], lwd=1.5)
segments(55, eac_s$s_55_lts[2], 60, eac_s$s_60_lts[2], col=brew[2], lwd=1.5)
#LHMPO
lines(5, eac_s$s_05_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(10, eac_s$s_10_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(15, eac_s$s_15_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(20, eac_s$s_20_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(25, eac_s$s_25_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(30, eac_s$s_30_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(35, eac_s$s_35_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(40, eac_s$s_40_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(45, eac_s$s_45_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(50, eac_s$s_50_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(55, eac_s$s_55_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(60, eac_s$s_60_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
segments(0, 0, 5, eac_s$s_05_lts[3], col=brew[3], lwd=1.5)
segments(5,  eac_s$s_05_lts[3], 10, eac_s$s_10_lts[3], col=brew[3], lwd=1.5)
segments(10, eac_s$s_10_lts[3], 15, eac_s$s_15_lts[3], col=brew[3], lwd=1.5)
segments(15, eac_s$s_15_lts[3], 20, eac_s$s_20_lts[3], col=brew[3], lwd=1.5)
segments(20, eac_s$s_20_lts[3], 25, eac_s$s_25_lts[3], col=brew[3], lwd=1.5)
segments(25, eac_s$s_25_lts[3], 30, eac_s$s_30_lts[3], col=brew[3], lwd=1.5)
segments(30, eac_s$s_30_lts[3], 35, eac_s$s_35_lts[3], col=brew[3], lwd=1.5)
segments(35, eac_s$s_35_lts[3], 40, eac_s$s_40_lts[3], col=brew[3], lwd=1.5)
segments(40, eac_s$s_40_lts[3], 45, eac_s$s_45_lts[3], col=brew[3], lwd=1.5)
segments(45, eac_s$s_45_lts[3], 50, eac_s$s_50_lts[3], col=brew[3], lwd=1.5)
segments(50, eac_s$s_50_lts[3], 55, eac_s$s_55_lts[3], col=brew[3], lwd=1.5)
segments(55, eac_s$s_55_lts[3], 60, eac_s$s_60_lts[3], col=brew[3], lwd=1.5)
#MAG
lines(5, eac_s$s_05_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(10, eac_s$s_10_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(15, eac_s$s_15_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(20, eac_s$s_20_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(25, eac_s$s_25_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(30, eac_s$s_30_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(35, eac_s$s_35_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(40, eac_s$s_40_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(45, eac_s$s_45_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(50, eac_s$s_50_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(55, eac_s$s_55_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(60, eac_s$s_60_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
segments(0, 0, 5, eac_s$s_05_lts[4], col=brew[4], lwd=1.5)
segments(5,  eac_s$s_05_lts[4], 10, eac_s$s_10_lts[4], col=brew[4], lwd=1.5)
segments(10, eac_s$s_10_lts[4], 15, eac_s$s_15_lts[4], col=brew[4], lwd=1.5)
segments(15, eac_s$s_15_lts[4], 20, eac_s$s_20_lts[4], col=brew[4], lwd=1.5)
segments(20, eac_s$s_20_lts[4], 25, eac_s$s_25_lts[4], col=brew[4], lwd=1.5)
segments(25, eac_s$s_25_lts[4], 30, eac_s$s_30_lts[4], col=brew[4], lwd=1.5)
segments(30, eac_s$s_30_lts[4], 35, eac_s$s_35_lts[4], col=brew[4], lwd=1.5)
segments(35, eac_s$s_35_lts[4], 40, eac_s$s_40_lts[4], col=brew[4], lwd=1.5)
segments(40, eac_s$s_40_lts[4], 45, eac_s$s_45_lts[4], col=brew[4], lwd=1.5)
segments(45, eac_s$s_45_lts[4], 50, eac_s$s_50_lts[4], col=brew[4], lwd=1.5)
segments(50, eac_s$s_50_lts[4], 55, eac_s$s_55_lts[4], col=brew[4], lwd=1.5)
segments(55, eac_s$s_55_lts[4], 60, eac_s$s_60_lts[4], col=brew[4], lwd=1.5)
#PAG
lines(5, eac_s$s_05_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(10, eac_s$s_10_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(15, eac_s$s_15_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(20, eac_s$s_20_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(25, eac_s$s_25_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(30, eac_s$s_30_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(35, eac_s$s_35_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(40, eac_s$s_40_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(45, eac_s$s_45_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(50, eac_s$s_50_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(55, eac_s$s_55_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(60, eac_s$s_60_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
segments(0, 0, 5, eac_s$s_05_lts[5], col=brew[5], lwd=1.5)
segments(5,  eac_s$s_05_lts[5], 10, eac_s$s_10_lts[5], col=brew[5], lwd=1.5)
segments(10, eac_s$s_10_lts[5], 15, eac_s$s_15_lts[5], col=brew[5], lwd=1.5)
segments(15, eac_s$s_15_lts[5], 20, eac_s$s_20_lts[5], col=brew[5], lwd=1.5)
segments(20, eac_s$s_20_lts[5], 25, eac_s$s_25_lts[5], col=brew[5], lwd=1.5)
segments(25, eac_s$s_25_lts[5], 30, eac_s$s_30_lts[5], col=brew[5], lwd=1.5)
segments(30, eac_s$s_30_lts[5], 35, eac_s$s_35_lts[5], col=brew[5], lwd=1.5)
segments(35, eac_s$s_35_lts[5], 40, eac_s$s_40_lts[5], col=brew[5], lwd=1.5)
segments(40, eac_s$s_40_lts[5], 45, eac_s$s_45_lts[5], col=brew[5], lwd=1.5)
segments(45, eac_s$s_45_lts[5], 50, eac_s$s_50_lts[5], col=brew[5], lwd=1.5)
segments(50, eac_s$s_50_lts[5], 55, eac_s$s_55_lts[5], col=brew[5], lwd=1.5)
segments(55, eac_s$s_55_lts[5], 60, eac_s$s_60_lts[5], col=brew[5], lwd=1.5)
#SCMPO
lines(5, eac_s$s_05_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(10, eac_s$s_10_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(15, eac_s$s_15_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(20, eac_s$s_20_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(25, eac_s$s_25_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(30, eac_s$s_30_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(35, eac_s$s_35_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(40, eac_s$s_40_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(45, eac_s$s_45_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(50, eac_s$s_50_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(55, eac_s$s_55_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(60, eac_s$s_60_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
segments(0, 0, 5, eac_s$s_05_lts[6], col=brew[6], lwd=1.5)
segments(5,  eac_s$s_05_lts[6], 10, eac_s$s_10_lts[6], col=brew[6], lwd=1.5)
segments(10, eac_s$s_10_lts[6], 15, eac_s$s_15_lts[6], col=brew[6], lwd=1.5)
segments(15, eac_s$s_15_lts[6], 20, eac_s$s_20_lts[6], col=brew[6], lwd=1.5)
segments(20, eac_s$s_20_lts[6], 25, eac_s$s_25_lts[6], col=brew[6], lwd=1.5)
segments(25, eac_s$s_25_lts[6], 30, eac_s$s_30_lts[6], col=brew[6], lwd=1.5)
segments(30, eac_s$s_30_lts[6], 35, eac_s$s_35_lts[6], col=brew[6], lwd=1.5)
segments(35, eac_s$s_35_lts[6], 40, eac_s$s_40_lts[6], col=brew[6], lwd=1.5)
segments(40, eac_s$s_40_lts[6], 45, eac_s$s_45_lts[6], col=brew[6], lwd=1.5)
segments(45, eac_s$s_45_lts[6], 50, eac_s$s_50_lts[6], col=brew[6], lwd=1.5)
segments(50, eac_s$s_50_lts[6], 55, eac_s$s_55_lts[6], col=brew[6], lwd=1.5)
segments(55, eac_s$s_55_lts[6], 60, eac_s$s_60_lts[6], col=brew[6], lwd=1.5)
#SVMPO
lines(5, eac_s$s_05_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(10, eac_s$s_10_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(15, eac_s$s_15_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(20, eac_s$s_20_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(25, eac_s$s_25_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(30, eac_s$s_30_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(35, eac_s$s_35_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(40, eac_s$s_40_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(45, eac_s$s_45_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(50, eac_s$s_50_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(55, eac_s$s_55_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(60, eac_s$s_60_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
segments(0, 0, 5, eac_s$s_05_lts[7], col=brew[7], lwd=1.5)
segments(5,  eac_s$s_05_lts[7], 10, eac_s$s_10_lts[7], col=brew[7], lwd=1.5)
segments(10, eac_s$s_10_lts[7], 15, eac_s$s_15_lts[7], col=brew[7], lwd=1.5)
segments(15, eac_s$s_15_lts[7], 20, eac_s$s_20_lts[7], col=brew[7], lwd=1.5)
segments(20, eac_s$s_20_lts[7], 25, eac_s$s_25_lts[7], col=brew[7], lwd=1.5)
segments(25, eac_s$s_25_lts[7], 30, eac_s$s_30_lts[7], col=brew[7], lwd=1.5)
segments(30, eac_s$s_30_lts[7], 35, eac_s$s_35_lts[7], col=brew[7], lwd=1.5)
segments(35, eac_s$s_35_lts[7], 40, eac_s$s_40_lts[7], col=brew[7], lwd=1.5)
segments(40, eac_s$s_40_lts[7], 45, eac_s$s_45_lts[7], col=brew[7], lwd=1.5)
segments(45, eac_s$s_45_lts[7], 50, eac_s$s_50_lts[7], col=brew[7], lwd=1.5)
segments(50, eac_s$s_50_lts[7], 55, eac_s$s_55_lts[7], col=brew[7], lwd=1.5)
segments(55, eac_s$s_55_lts[7], 60, eac_s$s_60_lts[7], col=brew[7], lwd=1.5)
#YMPO
lines(5, eac_s$s_05_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(10, eac_s$s_10_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(15, eac_s$s_15_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(20, eac_s$s_20_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(25, eac_s$s_25_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(30, eac_s$s_30_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(35, eac_s$s_35_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(40, eac_s$s_40_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(45, eac_s$s_45_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(50, eac_s$s_50_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(55, eac_s$s_55_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(60, eac_s$s_60_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
segments(0, 0, 5, eac_s$s_05_lts[8], col=brew[8], lwd=1.5)
segments(5,  eac_s$s_05_lts[8], 10, eac_s$s_10_lts[8], col=brew[8], lwd=1.5)
segments(10, eac_s$s_10_lts[8], 15, eac_s$s_15_lts[8], col=brew[8], lwd=1.5)
segments(15, eac_s$s_15_lts[8], 20, eac_s$s_20_lts[8], col=brew[8], lwd=1.5)
segments(20, eac_s$s_20_lts[8], 25, eac_s$s_25_lts[8], col=brew[8], lwd=1.5)
segments(25, eac_s$s_25_lts[8], 30, eac_s$s_30_lts[8], col=brew[8], lwd=1.5)
segments(30, eac_s$s_30_lts[8], 35, eac_s$s_35_lts[8], col=brew[8], lwd=1.5)
segments(35, eac_s$s_35_lts[8], 40, eac_s$s_40_lts[8], col=brew[8], lwd=1.5)
segments(40, eac_s$s_40_lts[8], 45, eac_s$s_45_lts[8], col=brew[8], lwd=1.5)
segments(45, eac_s$s_45_lts[8], 50, eac_s$s_50_lts[8], col=brew[8], lwd=1.5)
segments(50, eac_s$s_50_lts[8], 55, eac_s$s_55_lts[8], col=brew[8], lwd=1.5)
segments(55, eac_s$s_55_lts[8], 60, eac_s$s_60_lts[8], col=brew[8], lwd=1.5)
# Create x-axis
axis(side=1, at=c(0,5,10,15,20,25,30,35,40,45,50,55,60), tick=TRUE, cex=1.2)
# Create legend
text(0, 25, "Legend: EAC Cyclist", adj=0, cex=0.8, font=2, col="black")
points(0, 24, col=brew[1], cex=0.8, pch=20)
segments(0, 24, 5, 24, col=brew[1], lty=1, lwd=1.5)
text(7, 24, "CYMPO", adj=0, cex=0.8, col="black")
points(0, 23.5, col=brew[3], cex=0.8, pch=20)
segments(0, 23.5, 5, 23.5, col=brew[3], lty=1, lwd=1.5)
text(7, 23.5, "LHMPO", adj=0, cex=0.8, col="black")
points(0, 23, col=brew[4], cex=0.8, pch=20)
segments(0, 23, 5, 23, col=brew[4], lty=1, lwd=1.5)
text(7, 23, "MAG", adj=0, cex=0.8, col="black")
points(0, 22.5, col=brew[2], cex=0.8, pch=20)
segments(0, 22.5, 5, 22.5, col=brew[2], lty=1, lwd=1.5)
text(7, 22.5, "MetroPlan", adj=0, cex=0.8, col="black")
points(0, 22, col=brew[5], cex=0.8, pch=20)
segments(0, 22, 5, 22, col=brew[5], lty=1, lwd=1.5)
text(7, 22, "PAG", adj=0, cex=0.8, col="black")
points(0, 21.5, col=brew[6], cex=0.8, pch=20)
segments(0, 21.5, 5, 21.5, col=brew[6], lty=1, lwd=1.5)
text(7, 21.5, "SCMPO", adj=0, cex=0.8, col="black")
points(0, 21, col=brew[7], cex=0.8, pch=20)
segments(0, 21, 5, 21, col=brew[7], lty=1, lwd=1.5)
text(7, 21, "SVMPO", adj=0, cex=0.8, col="black")
points(0, 20.5, col=brew[8], cex=0.8, pch=20)
segments(0, 20.5, 5, 20.5, col=brew[8], lty=1, lwd=1.5)
text(7, 20.5, "YMPO", adj=0, cex=0.8, col="black")
rect(-1, 20, 14, 25.5, border="black", col=NA)
graphics.off()

### SAF Cyclist Type

# Graphic set-up
png(file=paste(dat_dir2, "plot_saf-schools.png", sep="/_analysis/_plots/"), bg="white", width=9, height=9, units='in', res=750)
plot(c(0,60), c(0,25), type="n", xlab="Travel time (minutes)", ylab="School accessibility", xaxt="n", cex=1.2)
# Plot points and lines
# CYMPO
lines(5, saf_s$s_05_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(10, saf_s$s_10_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(15, saf_s$s_15_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(20, saf_s$s_20_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(25, saf_s$s_25_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(30, saf_s$s_30_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(35, saf_s$s_35_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(40, saf_s$s_40_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(45, saf_s$s_45_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(50, saf_s$s_50_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(55, saf_s$s_55_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
lines(60, saf_s$s_60_lts[1], type="b", lwd=1.5, lty=0, col=brew[1], pch=20)
segments(0, 0, 5, saf_s$s_05_lts[1], col=brew[1], lwd=1.5)
segments(5,  saf_s$s_05_lts[1], 10, saf_s$s_10_lts[1], col=brew[1], lwd=1.5)
segments(10, saf_s$s_10_lts[1], 15, saf_s$s_15_lts[1], col=brew[1], lwd=1.5)
segments(15, saf_s$s_15_lts[1], 20, saf_s$s_20_lts[1], col=brew[1], lwd=1.5)
segments(20, saf_s$s_20_lts[1], 25, saf_s$s_25_lts[1], col=brew[1], lwd=1.5)
segments(25, saf_s$s_25_lts[1], 30, saf_s$s_30_lts[1], col=brew[1], lwd=1.5)
segments(30, saf_s$s_30_lts[1], 35, saf_s$s_35_lts[1], col=brew[1], lwd=1.5)
segments(35, saf_s$s_35_lts[1], 40, saf_s$s_40_lts[1], col=brew[1], lwd=1.5)
segments(40, saf_s$s_40_lts[1], 45, saf_s$s_45_lts[1], col=brew[1], lwd=1.5)
segments(45, saf_s$s_45_lts[1], 50, saf_s$s_50_lts[1], col=brew[1], lwd=1.5)
segments(50, saf_s$s_50_lts[1], 55, saf_s$s_55_lts[1], col=brew[1], lwd=1.5)
segments(55, saf_s$s_55_lts[1], 60, saf_s$s_60_lts[1], col=brew[1], lwd=1.5)
#FMPO
lines(5, saf_s$s_05_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(10, saf_s$s_10_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(15, saf_s$s_15_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(20, saf_s$s_20_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(25, saf_s$s_25_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(30, saf_s$s_30_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(35, saf_s$s_35_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(40, saf_s$s_40_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(45, saf_s$s_45_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(50, saf_s$s_50_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(55, saf_s$s_55_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
lines(60, saf_s$s_60_lts[2], type="b", lwd=1.5, lty=0, col=brew[2], pch=20)
segments(0, 0, 5, saf_s$s_05_lts[2], col=brew[2], lwd=1.5)
segments(5,  saf_s$s_05_lts[2], 10, saf_s$s_10_lts[2], col=brew[2], lwd=1.5)
segments(10, saf_s$s_10_lts[2], 15, saf_s$s_15_lts[2], col=brew[2], lwd=1.5)
segments(15, saf_s$s_15_lts[2], 20, saf_s$s_20_lts[2], col=brew[2], lwd=1.5)
segments(20, saf_s$s_20_lts[2], 25, saf_s$s_25_lts[2], col=brew[2], lwd=1.5)
segments(25, saf_s$s_25_lts[2], 30, saf_s$s_30_lts[2], col=brew[2], lwd=1.5)
segments(30, saf_s$s_30_lts[2], 35, saf_s$s_35_lts[2], col=brew[2], lwd=1.5)
segments(35, saf_s$s_35_lts[2], 40, saf_s$s_40_lts[2], col=brew[2], lwd=1.5)
segments(40, saf_s$s_40_lts[2], 45, saf_s$s_45_lts[2], col=brew[2], lwd=1.5)
segments(45, saf_s$s_45_lts[2], 50, saf_s$s_50_lts[2], col=brew[2], lwd=1.5)
segments(50, saf_s$s_50_lts[2], 55, saf_s$s_55_lts[2], col=brew[2], lwd=1.5)
segments(55, saf_s$s_55_lts[2], 60, saf_s$s_60_lts[2], col=brew[2], lwd=1.5)
#LHMPO
lines(5, saf_s$s_05_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(10, saf_s$s_10_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(15, saf_s$s_15_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(20, saf_s$s_20_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(25, saf_s$s_25_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(30, saf_s$s_30_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(35, saf_s$s_35_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(40, saf_s$s_40_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(45, saf_s$s_45_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(50, saf_s$s_50_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(55, saf_s$s_55_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
lines(60, saf_s$s_60_lts[3], type="b", lwd=1.5, lty=0, col=brew[3], pch=20)
segments(0, 0, 5, saf_s$s_05_lts[3], col=brew[3], lwd=1.5)
segments(5,  saf_s$s_05_lts[3], 10, saf_s$s_10_lts[3], col=brew[3], lwd=1.5)
segments(10, saf_s$s_10_lts[3], 15, saf_s$s_15_lts[3], col=brew[3], lwd=1.5)
segments(15, saf_s$s_15_lts[3], 20, saf_s$s_20_lts[3], col=brew[3], lwd=1.5)
segments(20, saf_s$s_20_lts[3], 25, saf_s$s_25_lts[3], col=brew[3], lwd=1.5)
segments(25, saf_s$s_25_lts[3], 30, saf_s$s_30_lts[3], col=brew[3], lwd=1.5)
segments(30, saf_s$s_30_lts[3], 35, saf_s$s_35_lts[3], col=brew[3], lwd=1.5)
segments(35, saf_s$s_35_lts[3], 40, saf_s$s_40_lts[3], col=brew[3], lwd=1.5)
segments(40, saf_s$s_40_lts[3], 45, saf_s$s_45_lts[3], col=brew[3], lwd=1.5)
segments(45, saf_s$s_45_lts[3], 50, saf_s$s_50_lts[3], col=brew[3], lwd=1.5)
segments(50, saf_s$s_50_lts[3], 55, saf_s$s_55_lts[3], col=brew[3], lwd=1.5)
segments(55, saf_s$s_55_lts[3], 60, saf_s$s_60_lts[3], col=brew[3], lwd=1.5)
#MAG
lines(5, saf_s$s_05_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(10, saf_s$s_10_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(15, saf_s$s_15_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(20, saf_s$s_20_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(25, saf_s$s_25_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(30, saf_s$s_30_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(35, saf_s$s_35_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(40, saf_s$s_40_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(45, saf_s$s_45_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(50, saf_s$s_50_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(55, saf_s$s_55_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
lines(60, saf_s$s_60_lts[4], type="b", lwd=1.5, lty=0, col=brew[4], pch=20)
segments(0, 0, 5, saf_s$s_05_lts[4], col=brew[4], lwd=1.5)
segments(5,  saf_s$s_05_lts[4], 10, saf_s$s_10_lts[4], col=brew[4], lwd=1.5)
segments(10, saf_s$s_10_lts[4], 15, saf_s$s_15_lts[4], col=brew[4], lwd=1.5)
segments(15, saf_s$s_15_lts[4], 20, saf_s$s_20_lts[4], col=brew[4], lwd=1.5)
segments(20, saf_s$s_20_lts[4], 25, saf_s$s_25_lts[4], col=brew[4], lwd=1.5)
segments(25, saf_s$s_25_lts[4], 30, saf_s$s_30_lts[4], col=brew[4], lwd=1.5)
segments(30, saf_s$s_30_lts[4], 35, saf_s$s_35_lts[4], col=brew[4], lwd=1.5)
segments(35, saf_s$s_35_lts[4], 40, saf_s$s_40_lts[4], col=brew[4], lwd=1.5)
segments(40, saf_s$s_40_lts[4], 45, saf_s$s_45_lts[4], col=brew[4], lwd=1.5)
segments(45, saf_s$s_45_lts[4], 50, saf_s$s_50_lts[4], col=brew[4], lwd=1.5)
segments(50, saf_s$s_50_lts[4], 55, saf_s$s_55_lts[4], col=brew[4], lwd=1.5)
segments(55, saf_s$s_55_lts[4], 60, saf_s$s_60_lts[4], col=brew[4], lwd=1.5)
#PAG
lines(5, saf_s$s_05_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(10, saf_s$s_10_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(15, saf_s$s_15_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(20, saf_s$s_20_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(25, saf_s$s_25_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(30, saf_s$s_30_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(35, saf_s$s_35_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(40, saf_s$s_40_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(45, saf_s$s_45_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(50, saf_s$s_50_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(55, saf_s$s_55_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
lines(60, saf_s$s_60_lts[5], type="b", lwd=1.5, lty=0, col=brew[5], pch=20)
segments(0, 0, 5, saf_s$s_05_lts[5], col=brew[5], lwd=1.5)
segments(5,  saf_s$s_05_lts[5], 10, saf_s$s_10_lts[5], col=brew[5], lwd=1.5)
segments(10, saf_s$s_10_lts[5], 15, saf_s$s_15_lts[5], col=brew[5], lwd=1.5)
segments(15, saf_s$s_15_lts[5], 20, saf_s$s_20_lts[5], col=brew[5], lwd=1.5)
segments(20, saf_s$s_20_lts[5], 25, saf_s$s_25_lts[5], col=brew[5], lwd=1.5)
segments(25, saf_s$s_25_lts[5], 30, saf_s$s_30_lts[5], col=brew[5], lwd=1.5)
segments(30, saf_s$s_30_lts[5], 35, saf_s$s_35_lts[5], col=brew[5], lwd=1.5)
segments(35, saf_s$s_35_lts[5], 40, saf_s$s_40_lts[5], col=brew[5], lwd=1.5)
segments(40, saf_s$s_40_lts[5], 45, saf_s$s_45_lts[5], col=brew[5], lwd=1.5)
segments(45, saf_s$s_45_lts[5], 50, saf_s$s_50_lts[5], col=brew[5], lwd=1.5)
segments(50, saf_s$s_50_lts[5], 55, saf_s$s_55_lts[5], col=brew[5], lwd=1.5)
segments(55, saf_s$s_55_lts[5], 60, saf_s$s_60_lts[5], col=brew[5], lwd=1.5)
#SCMPO
lines(5, saf_s$s_05_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(10, saf_s$s_10_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(15, saf_s$s_15_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(20, saf_s$s_20_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(25, saf_s$s_25_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(30, saf_s$s_30_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(35, saf_s$s_35_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(40, saf_s$s_40_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(45, saf_s$s_45_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(50, saf_s$s_50_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(55, saf_s$s_55_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
lines(60, saf_s$s_60_lts[6], type="b", lwd=1.5, lty=0, col=brew[6], pch=20)
segments(0, 0, 5, saf_s$s_05_lts[6], col=brew[6], lwd=1.5)
segments(5,  saf_s$s_05_lts[6], 10, saf_s$s_10_lts[6], col=brew[6], lwd=1.5)
segments(10, saf_s$s_10_lts[6], 15, saf_s$s_15_lts[6], col=brew[6], lwd=1.5)
segments(15, saf_s$s_15_lts[6], 20, saf_s$s_20_lts[6], col=brew[6], lwd=1.5)
segments(20, saf_s$s_20_lts[6], 25, saf_s$s_25_lts[6], col=brew[6], lwd=1.5)
segments(25, saf_s$s_25_lts[6], 30, saf_s$s_30_lts[6], col=brew[6], lwd=1.5)
segments(30, saf_s$s_30_lts[6], 35, saf_s$s_35_lts[6], col=brew[6], lwd=1.5)
segments(35, saf_s$s_35_lts[6], 40, saf_s$s_40_lts[6], col=brew[6], lwd=1.5)
segments(40, saf_s$s_40_lts[6], 45, saf_s$s_45_lts[6], col=brew[6], lwd=1.5)
segments(45, saf_s$s_45_lts[6], 50, saf_s$s_50_lts[6], col=brew[6], lwd=1.5)
segments(50, saf_s$s_50_lts[6], 55, saf_s$s_55_lts[6], col=brew[6], lwd=1.5)
segments(55, saf_s$s_55_lts[6], 60, saf_s$s_60_lts[6], col=brew[6], lwd=1.5)
#SVMPO
lines(5, saf_s$s_05_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(10, saf_s$s_10_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(15, saf_s$s_15_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(20, saf_s$s_20_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(25, saf_s$s_25_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(30, saf_s$s_30_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(35, saf_s$s_35_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(40, saf_s$s_40_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(45, saf_s$s_45_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(50, saf_s$s_50_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(55, saf_s$s_55_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
lines(60, saf_s$s_60_lts[7], type="b", lwd=1.5, lty=0, col=brew[7], pch=20)
segments(0, 0, 5, saf_s$s_05_lts[7], col=brew[7], lwd=1.5)
segments(5,  saf_s$s_05_lts[7], 10, saf_s$s_10_lts[7], col=brew[7], lwd=1.5)
segments(10, saf_s$s_10_lts[7], 15, saf_s$s_15_lts[7], col=brew[7], lwd=1.5)
segments(15, saf_s$s_15_lts[7], 20, saf_s$s_20_lts[7], col=brew[7], lwd=1.5)
segments(20, saf_s$s_20_lts[7], 25, saf_s$s_25_lts[7], col=brew[7], lwd=1.5)
segments(25, saf_s$s_25_lts[7], 30, saf_s$s_30_lts[7], col=brew[7], lwd=1.5)
segments(30, saf_s$s_30_lts[7], 35, saf_s$s_35_lts[7], col=brew[7], lwd=1.5)
segments(35, saf_s$s_35_lts[7], 40, saf_s$s_40_lts[7], col=brew[7], lwd=1.5)
segments(40, saf_s$s_40_lts[7], 45, saf_s$s_45_lts[7], col=brew[7], lwd=1.5)
segments(45, saf_s$s_45_lts[7], 50, saf_s$s_50_lts[7], col=brew[7], lwd=1.5)
segments(50, saf_s$s_50_lts[7], 55, saf_s$s_55_lts[7], col=brew[7], lwd=1.5)
segments(55, saf_s$s_55_lts[7], 60, saf_s$s_60_lts[7], col=brew[7], lwd=1.5)
#YMPO
lines(5, saf_s$s_05_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(10, saf_s$s_10_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(15, saf_s$s_15_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(20, saf_s$s_20_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(25, saf_s$s_25_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(30, saf_s$s_30_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(35, saf_s$s_35_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(40, saf_s$s_40_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(45, saf_s$s_45_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(50, saf_s$s_50_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(55, saf_s$s_55_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
lines(60, saf_s$s_60_lts[8], type="b", lwd=1.5, lty=0, col=brew[8], pch=20)
segments(0, 0, 5, saf_s$s_05_lts[8], col=brew[8], lwd=1.5)
segments(5,  saf_s$s_05_lts[8], 10, saf_s$s_10_lts[8], col=brew[8], lwd=1.5)
segments(10, saf_s$s_10_lts[8], 15, saf_s$s_15_lts[8], col=brew[8], lwd=1.5)
segments(15, saf_s$s_15_lts[8], 20, saf_s$s_20_lts[8], col=brew[8], lwd=1.5)
segments(20, saf_s$s_20_lts[8], 25, saf_s$s_25_lts[8], col=brew[8], lwd=1.5)
segments(25, saf_s$s_25_lts[8], 30, saf_s$s_30_lts[8], col=brew[8], lwd=1.5)
segments(30, saf_s$s_30_lts[8], 35, saf_s$s_35_lts[8], col=brew[8], lwd=1.5)
segments(35, saf_s$s_35_lts[8], 40, saf_s$s_40_lts[8], col=brew[8], lwd=1.5)
segments(40, saf_s$s_40_lts[8], 45, saf_s$s_45_lts[8], col=brew[8], lwd=1.5)
segments(45, saf_s$s_45_lts[8], 50, saf_s$s_50_lts[8], col=brew[8], lwd=1.5)
segments(50, saf_s$s_50_lts[8], 55, saf_s$s_55_lts[8], col=brew[8], lwd=1.5)
segments(55, saf_s$s_55_lts[8], 60, saf_s$s_60_lts[8], col=brew[8], lwd=1.5)
# Create x-axis
axis(side=1, at=c(0,5,10,15,20,25,30,35,40,45,50,55,60), tick=TRUE, cex=1.2)
# Create legend
text(0, 25, "Legend: SAF Cyclist", adj=0, cex=0.8, font=2, col="black")
points(0, 24, col=brew[1], cex=0.8, pch=20)
segments(0, 24, 5, 24, col=brew[1], lty=1, lwd=1.5)
text(7, 24, "CYMPO", adj=0, cex=0.8, col="black")
points(0, 23.5, col=brew[3], cex=0.8, pch=20)
segments(0, 23.5, 5, 23.5, col=brew[3], lty=1, lwd=1.5)
text(7, 23.5, "LHMPO", adj=0, cex=0.8, col="black")
points(0, 23, col=brew[4], cex=0.8, pch=20)
segments(0, 23, 5, 23, col=brew[4], lty=1, lwd=1.5)
text(7, 23, "MAG", adj=0, cex=0.8, col="black")
points(0, 22.5, col=brew[2], cex=0.8, pch=20)
segments(0, 22.5, 5, 22.5, col=brew[2], lty=1, lwd=1.5)
text(7, 22.5, "MetroPlan", adj=0, cex=0.8, col="black")
points(0, 22, col=brew[5], cex=0.8, pch=20)
segments(0, 22, 5, 22, col=brew[5], lty=1, lwd=1.5)
text(7, 22, "PAG", adj=0, cex=0.8, col="black")
points(0, 21.5, col=brew[6], cex=0.8, pch=20)
segments(0, 21.5, 5, 21.5, col=brew[6], lty=1, lwd=1.5)
text(7, 21.5, "SCMPO", adj=0, cex=0.8, col="black")
points(0, 21, col=brew[7], cex=0.8, pch=20)
segments(0, 21, 5, 21, col=brew[7], lty=1, lwd=1.5)
text(7, 21, "SVMPO", adj=0, cex=0.8, col="black")
points(0, 20.5, col=brew[8], cex=0.8, pch=20)
segments(0, 20.5, 5, 20.5, col=brew[8], lty=1, lwd=1.5)
text(7, 20.5, "YMPO", adj=0, cex=0.8, col="black")
rect(-1, 20, 14, 25.5, border="black", col=NA)
graphics.off()

### NOTE: Create three figures (each destination, all cyclist types) with produced plots

####################
### END OF SCRIPT
